Generalized dynamical equation and low energy nucleon dynamics 
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Low energy nucleon dynamics is investigated by using the generalized dynamical equation derived 
in [J. Phys. A v. 32, 5657 (1999)]. This equation extends quantum dynamics to describe the time 
evolution in the case of nonlocal-in-time interactions. We show that the use of the generalized 
dynamical equation, which allows one to take into account that effective actions in quantum field 
theory, arising after integrating out certain degrees of freedom, are generally nonlocal, provides a 
new way to formulate the effective theory of nuclear forces. In this way we construct the most 
general possible operator describing the short-range part of the NN interaction in the ^5*0 channel. 
In particular, this operator can be incorporated into the chiral potential model. This is shown to 
provide an extension of the standard chiral approach. Being equivalent to the chiral potential model 
for some values of the constants characterizing the short-range of the NN interaction this approach 
allows the solutions that cannot be reproduced by any potential model. In this way we find that 
for certain values of the constants inconsistent with potentials models the effect of the short-range 
forces on the half off-shell 27V T-matrix with high off-shell momenta grows in magnitude rapidly as 
the kinetic energy of outgoing protons decreases at laboratory energies below 15 MeV. We show 
that this effect may be the origin of the existing discrepancy between the theory and experiment in 
describing the proton-proton bremsstrahlung. 



I. INTRODUCTION 



Ideas from the foundations of quantum mechanics are 
being applied now to many branches of physics. In the 
quantum mechanics of particles interacting through the 
Coulomb potential one deals with a well-defined inter- 
action Hamiltonian and the Schrodinger equation gov- 
erning the dynamics of the theory. This theory is per- 
fectly consistent and provides an excellent description of 
atomic phenomena at low energies. It is natural to ex- 
pect that low energy nuclear physics can be described 
in the same way. However, one has not yet constructed 
a fundamental nucleon-nucleon (NN) potential. Nowa- 
days there exist phenomenological NN potentials which 
successfully describe scattering data to high precision, 
but they do not emerge from QCD and contain ad hoc 
form factors. A first attempt to systematically solve the 
problem of low energy nucleon dynamics and construct a 
bridge to QCD was made by Weinberg 1]. He suggested 
to derive a NN potential in time-ordered chiral perturba- 
tion theory (ChPT). However, such a potential is singu- 
lar and the Schrodinger (Lippmann-Schwinger) equation 
does not make sense without regularization and renor- 
malization. This means that in the effective field theory 
(EFT) of nuclear forces, which following the pioneering 
work of Weinberg has become very popular in nuclear 
physics (for a review, see Ref. 0), the Schrodinger equa- 
tion is not valid. On the other hand, the whole formalism 
of fields and particles can be considered as an inevitable 
consequence of quantum mechanics, Lorentz invariance, 
and the cluster decomposition principle 3j. Thus in the 
nonrelativistic limit QCD must reproduce low energy nu- 
cleon physics consistent with the basic principles of quan- 
tum mechanics. However, as it follows from the Weinberg 
analysis, QCD leads through ChPT to the low energy 
theory in which the Schrodinger equation is not valid. 



This means that either there is something wrong with 
QCD and ChPT or the Schrodinger equation is not the 
basic dynamical equation of quantum theory. Meanwhile, 
in Ref. Q it has been shown that the Schrodinger equa- 
tion is not the most general equation consistent with the 
current concepts of quantum physics, and a more general 
equation of motion has been derived as a consequence 
of the basic postulates of the Feynman and canonical 
approaches to quantum theory. Being equivalent to the 
Schrodinger equation in the case of instantaneous inter- 
actions, this generalized dynamical equation permits the 
generalization to the case where the dynamics of a sys- 
tem is generated by a nonlocal-in-time interaction. The 
generalized quantum dynamics (GQD) developed in this 
way has proved an useful tool for solving various prob- 
lems in quantum theory Q ■ 

The GQD allows one to consider the problem of con- 
sistency of quantum mechanics with the low energy pre- 
dictions of QCD from a new point of view. From this 
viewpoint, in investigating the consequences of ChPT we 
must not restrict ourselves to the assumption that the 
NN interaction can be parametrized by a NN poten- 
tial, and low energy nucleon dynamics is governed by the 
Schrodinger equation. This dynamics may be governed 
by the generalized dynamical equation (GDE) with a 
nonlocal-in-time interaction operator when this equation 
is not equivalent to the Schrodinger equation, and hence 
the above divergence problems may be the cost of try- 
ing to describe low energy nucleon dynamics in terms of 
Hamiltonian formalism while this dynamics is really non- 
Hamiltonian. Another motivation for investigating the 
problem of the description of nucleon dynamics from the 
point of view of the formalism of the GQD is still existing 
discrepancy between the predictions of the modern NN 
potential models and experiment. For example, there is 
a significant discrepancy between theory and experiment 
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in describing the proton-proton (pp) bremsstrahlung Q 
that is an accurate test of NN interaction models. An 
important part of this discrepancy originates in a poor 
description of the NN interaction at low energies 
The formalism of the GQD that reproduces the potential 
models of the NN interaction in the particular case when 
this interaction is assumed to be instantaneous admits a 
more general class of such models, and only experiment 
such as the pp bremsstrahlung can discriminate between 
them. In this context the discrepancy between the pre- 
dictions of the potential models and experiment should 
mean that in describing nucleon dynamics one cannot 
ignore that the NN interaction is nonlocal in time. 

The aim of the present paper is twofold. From the 
one hand, we intend to show that the low energy predic- 
tions of QCD are consistent with the basic principles of 
quantum mechanics, but in this case a new insight into 
these principles provided by the formalism of the GQD is 
needed. On the other hand, it is our intention to demon- 
strate that this formalism allows one to formulate the 
effective theory of nuclear forces as an inevitable conse- 
quence of the basic principles of quantum mechanics and 
the symmetries of QCD. Being formulated in this way 
the theory of nuclear forces opens new possibilities for 
describing nucleon dynamics at low energies. In Sec. II 
we briefly consider the main features of the formalism 
of the GQD developed by one of the authors (R.G.) in 
Ref. 5| ■ We mainly focus on the physical meaning of the 
GDE which is a direct consequence of the principle of 
the superposition of the probability amplitudes and the 
requirement that the evolution operator is unitary. This 
equation plays a key role in the present work. 

In Sec. Ill we consider the pionless theory of nuclear 
forces. By focusing the attention on the ^5*0 channel, we 
show that the requirement that the two-nucleon (2iV) T- 
matrix satisfies the GDE and is consistent with the sym- 
metries of QCD allows one to construct this T-matrix by 
expanding it in powers of Q/A where Q being some low 
energy scale, and A is the scale at which the theory is 
expected to break down. In this way we can reproduce 
all results for the NN scattering amplitudes obtained 
in the standard EFT of nuclear forces. It is important 
that these results are reproduced starting with a well 
defined 2N T-matrix without resorting to the regular- 
ization and renormalization procedures. In Sec. IV we 
show that the use of the GDE allows one to formulate 
the effective theory as a perfectly consistent theory free 
from UV divergences. Being formulated in this way, the 
effective theory keeps all advantages of the traditional 
nuclear physics approach. At the same time, its advan- 
tage over the traditional approach is that it allows one 
to find constraints on the off-shell behavior of the 2N 
T-matrix that are placed by the symmetries of QCD. 

In Sec. V the proposed formalism is considered from 
the point of view of the Weinberg program for physics 
of the two-nucleon systems. The approach based on the 
derivation of a NN potential from ChPT was pioneered 
by Weinberg |^ and developed by Ordonez ^3 ^.nd van 



Kolck Recently the potential at fourth order of 

ChPT that reproduce the NN data with the same accu- 
racy as phenomenological high-precision potentials was 
developed by Entem and Machleidt T^. We show that 
our formulation of the effective theory of nuclear forces 
is a generalization of the above approach. Instead of 
the Schrodinger (LS) equation, it is suggested to use the 
GDE. In this way we construct the most general possible 
operator describing the short-range part of the NN inter- 
action. We show that this operator can be incorporated 
into the chiral approach. Being equivalent to the chiral 
potential model for some values of the constants charac- 
terizing the short-range part of the NN interaction this 
approach allows the solutions that cannot be reproduced 
by any potential model. In Sec. VI we show that the 
for certain values of the constants inconsistent with po- 
tentials models the effect of the short-range forces on the 
2A'^ T-matrix with high off-shell momenta grows in mag- 
nitude rapidly as the kinetic energy of outgoing protons 
decreases in the region below 15 MeV. We also show that 
the fact that the existing microscopic pp bremsstrahlung 
models cannot reproduce this effect may be the origin of 
the discrepancy between their predictions and the exper- 
imental data. We also show that our formalism allows 
one to modify the short-range parts of the existing high- 
precision potentials to improve their off-shell predictions 
keeping the phase shifts unchanged. New possibilities 
that the formalism of the GQD opens for solving three- 
body problem in nuclear physics are discussed in Sec. 
VII. 



II. GENERALIZED QUANTUM DYNAMICS 

The basic concept of the canonical formalism of quan- 
tum theory is that it can be formulated in terms of 
vectors of a Hilbert space and operators acting on this 
space. In this formalism the postulates that establish 
the connection between the vectors and operators and 
states of a quantum system and observables are used to- 
gether with the dynamical postulate according to which 
the time evolution of a quantum system is governed by 
the Schrodinger equation. In the Feynman formalism 
quantum theory is formulated in terms of probability 
amplitudes without resorting to vectors and operators 
acting on a Hilbert space. Feynman's theory starts with 
an analysis of the phenomenon of quantum interference. 
The results of this analysis which leads directly to the 
concept of the superposition of probability amplitudes 
are summarized by the following postulate 5j: 

The probability of an event is the absolute square of 
a complex number called the probability amplitude. The 
joint probability amplitude of a time-ordered sequence of 
events is the product of separate probability amplitudes 
of each of these events. The probability amplitude of an 
event which can happen in several different ways is a sum 
of the probability amplitudes for each of these ways. 

The Feynman formulation is based on the assumption 
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that the history of a quantum system can be represented 
by some path in space-time, and hence the probabihty 
amphtude of any event is a sum of the probabihty amph- 
tudes that a particle has a completely specified path in 
space-time. The contribution from a single path is postu- 
lated to be an exponential whose (imaginary) phase is the 
classical action (in units of h) for the path in question. 
This assumption is not as fundamental as the above prin- 
ciple of the superposition of probability amplitudes which 
follows directly from the analyzes of the phenomenon of 
quantum interference. This fact is emphasized, for ex- 
ample, in Feynman's book |l3l |. where a minimal set of 
physical principles which must be satisfied in any theory 
of fields and particles is analyzed. Feynman includes in 
this set neither the second postulate of his formalism nor 
the Schrodinger equation: The only quantum mechanical 
principle included in this set is the principle of the su- 
perposition of probability amplitudes. In Ref. I3| it has 
been shown that, instead of processes associated with a 
completely specified path in space-time, one can use pro- 
cesses associated with completely specified instants of the 
beginning and end of interaction in a quantum system as 
alternative ways in which any event can happen. As it 
turned out, employing this class of alternatives allows 
one to derive a dynamical equation from the principle 
of the superposition without making supplementary as- 
sumptions like the second postulate of Feynman's theory. 
By using this class of alternatives and the superposition 
principle, {ip\U{t,to)\ip) ^ being the probability amplitude 
of finding a quantum system in the state \ip) in a mea- 
surement at time t, if at time it was in the state \(p), 
can be represented in the form |j| 

{^p\u{t,toM^{i^\v)+ I dt2 r dh{i^\~s{t2Mmii) 

Jto Jta 

Here ('0|5'(<2, ii)|<(5) is the probability amplitude that, if 
at time ti the system was in the state \lp), then the in- 
teraction in the system will begin at time ti and end at 
time t2 and at this time the system will be in the state 
IV'). By using the operator formalism, one can represent 
amplitudes {'4>\U{t,t^))\Lp) by the matrix elements of the 
unitary evolution operator U {t, Iq) in the interaction pic- 
ture. The operator S{t2,ti) represents the contribution 
to the evolution operator from the process in which the 
interaction in the system begins at time ti and ends at 
time t2 ■ As has been shown in Ref. , for the evolution 
operator JQ) to be unitary for any t and to the operator 
5(^2,^1) must satisfy the equation 

it2-ti)S{t2,ti) 

= / dU dt3iU~t3)Sit2,U)Sit3,h). (2) 
Jtl Jtl 

A remarkable feature of this equation is that it works as 
a recursion relation and allows one to obtain the opera- 
tors S{t2,ti) for any ti and ^2, if S{t'2,t'i) corresponding 
to infinitesimal duration times r = ~ t'l of interaction 



are known. It is natural to assume that most of the con- 
tribution to the evolution operator in the limit t2 — > ii 
comes from the processes associated with the fundamen- 
tal interaction in the system under study. Denoting this 
contribution by Hint{t2,ii) we can write 

S{t2M) ^ H,^t{t2M)+o{T'), (3) 

where t = t2 — ti. The parameter e is determined by 
demanding that i?int (^2,^1) called the generalized inter- 
action operator must be so close to the solution of Eq. 

in the limit t2 — > ii that this equation has a unique 
solution having the behavior (j2Jl near the point ^2 = ^i- If 
Hintih, ti) is specified, Eq. Q allows one to find the op- 
erator 5(^2,^1), and hence the evolution operator. Thus 
Eq. (|2Jl can be regarded as an equation of motion for 
states of a quantum system. This generalized dynamical 
equation allows one to construct the evolution operator 
by using the contributions from fundamental processes 
as building blocks. In the case of Hamiltonian dynamics 
the fundamental interaction is instantaneous. The gener- 
alized interaction operator describing such an interaction 
is of the form 

HMit2,ti) = -2^(5(^2 - h)Hi{ti) (4) 

(the delta function 6(t2^ti) emphasizes that the interac- 
tion is instantaneous). In this case Eq. ^ is equivalent 
to the Schrodinger equation jj] with the operator Hj{t) 
being an interaction Hamiltonian (in the interaction pic- 
ture). At the same time, Eq. ||2J) permits the gener- 
alization to the case where the fundamental interaction 
in a quantum system is nonlocal in time, and hence the 
dynamics is non-Hamiltonian Q. 

Let us now briefly consider (for details see Ref. 3) the 
mathematical assumptions on which the formalism of the 
GQD is founded. As we already noted, in the canonical 
formalism the assumption that the evolution of a quan- 
tum system is governed by the Schrodinger equation is 
used as a dynamical postulate. On the other hand, since 
the operator V{t) = U{t, 0) is postulated to be unitary 
with the group properties 

V{t2,h)^Vit2)V{h), 1/(0) = 1, 

from Stone's theorem it follows that, if this one- 
parameter group is weekly continuous, i.e., 

(^-2 1^*2)1^-1),^, (^2|^(ii)|^i) (5) 

for any and \ip2) belonging to the Hilbert space, than 
it has a self-adjoint infinitesimal generator H: 

V{t) = exp{-iHt), idV{t)/dt = HV{t). 

Identifying H with the total Hamiltonian as usual, we get 
the Schrodinger equation. Thus the assumption that this 
equation governs the dynamics of a system is equivalent 
to the assumption that the evolution operator satisfies 
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the condition However, from the physical point of 
view, Eq. must not be satisfied for any states belong- 
ing to the Hilbert space. In fact, there are normalized 
vectors in the Hilbert space that represent the states for 
which energy of a system is infinite. From the point of 
view of the states with infinite energy |V'm) any time 
interval St is infinite, and hence the evolution operator 
U{St,0) = 1 + R{St,0) must be independent of 6t, i.e., 
must be constant 

(V','jc/(<5i,0)|v™) = (^^J^„) + (V4J^(^,0)|^,„). 

This means that for the evolution operator to be weekly 
continuous, {^[j^\R{oo,0)\'4'in) must be zero. However, 
in quantum field theory, for example, this is not the case 
because in this theory renormalization is necessary to ef- 
fectively cut off intermediate states with infinite energy. 
On the other hand, the states \tpin) are not physically 
realizable, and hence from the physical point of view it is 
enough to require that the condition jSJl must be satisfied 
only for the physically realizable states. The main idea 
of the formalism of the GQD is that in describing the 
dynamics of a system one need not to restrict oneself to 
the assumption that the evolution operator is necessarily 
weekly continuous, and therefore dynamics is necessarily 
governed by the Schrodinger equation. In Ref. j^] it has 
been shown that the current concepts of quantum physics 
generate the dynamical principle in which the above as- 
sumption does not play any role. In this way the GDE 
has been derived as the most general dynamical equation. 
This equation allows one to construct the evolution op- 
erator starting with the contributions from the processes 
with infinitesimal duration times of interaction that are 
described by the interaction operator i?mt (^2,^1), while 
the Schrodinger equation constructs the evolution oper- 
ator starting with the contributions from the processes 
with instantaneous interaction that are described by the 
interaction Hamiltonian. The GDE allows one to de- 
scribe, in a consistent way, the dynamics of a quantum 
system with a nonlocal in time interaction. This equa- 
tion is equivalent to the Schrodinger equation in the case 
when the interaction in a system is instantaneous. In 
the case of nonlocal-in-time interactions the dynamics 
is determined by the behavior of the interaction operator 
Hint{t2, ti) in the infinitesimal neighborhood of the point 
^2 = ii, i-e., such a nonlocality has no a definite scale. 
This scale should depend on the problem we consider: 
Some duration time of interaction may seem as infinites- 
imal from the point of view of low energy physics, but 
at the same time it may seem as finite from viewpoint of 
high energy physics. Thus the GDE allows one to take 
into account, in a consistent way, that every theory with 
which we deal is a low energy approximation to a more 
fundamental one and provides a bridge between them. 
This is may be important for describing the dynamics 
in the effective theory of nuclear forces. In this case the 
duration times of interaction much above the time scale 
of low energy nuclear physics but much below the scale 
at which high energy degrees of freedom come into play 



should be considered as "infinitesimal", and the operator 
describing the processes with such duration times of in- 
teraction can be used as an effective interaction operator. 
By using Eq. (QJ, for C/(i,to), we can write 



U{t,to) = 1 + — / dxexp[-i(z- iJo)i] (6) 
X (z - Ho)-^T{z){z - Ho)-' exp[i(z - Ho)tol 



where z =^ x + iy, y > 0, and Hq is the free Hamiltonian. 
The operator T(z) is defined by 



/•oo 

T{z) — i I dr exp{izT)T{T), 
Jo 



(7) 



where T(r) = exp(— ii/oi2)»S'(i2, ii) exp(iHoti), and r = 
t2 — ii . In terms of the T-matrix defined by Eq. Q the 
equation of motion ^ can be rewritten in the form Q 



d(V>2|r(^)|^i 

dz 



E 



{^2\T{z)\n){n\T{z)\^, 



(8) 



where n stands for the entire set of discrete and contin- 
uous variables that characterize the system in full, and 
\n) are the eigenvectors of Hq. As it follows from Eqs. 
(O and Q, the boundary condition on this equation is 
of the form 



{^P2\B{zM^) + o{\z\-^) 

DO 

{^2\B{z)\,Pi) + 0{hiz)}, (9) 



B{z) = t dTeM^^r)Htl{T), (10) 
^0 



where 



Hintih - h) ex.p{-iHot2)Hi„tit2,ti) exp{iHoti) 

is the interaction operator in the Schrodinger picture, 
and h{z) is an arbitrary function satisfying the condi- 



00, with P = e + 1. In 

is). 



tion h{z) — o [\z\ 

the case of the Hamiltonian dynamics when H^^'^ir) 
—2i6(T)Hj , with Hj being the interaction Hamiltonian 
in the Schrodinger picture, B{z) = Hj, and the boundary 
condition © takes the form 



T(z) ^ Hj. 

\z\ — > 00 



(11) 



Equation (jS)) with this boundary condition is equivalent 
to the Lippmann-Schwinger (LS) equation with the inter- 
action Hamiltonian Hj. By definition, the operator B{z) 
is the interaction operator in the energy representation. 
This operator must be so close to the relevant solution 
of Eq. in the limit \z\ ^00 that this differential 
equation has a unique solution having the asymptotic 
behavior (j5J. For this the operator B{z) must satisfy the 
condition 



d{^2\B{z)\^i 
dz 



^ {^2\B{z)\n){n\B{z)\iji) 



o{\z\ 



00. 



(12) 
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From Eq. JBJ it follows that the evolution operator in 
the Schrodinger picture can be represented in the form 



Us{t,0) = ^ y" dxexp{-izt)G{z), 



where z — x + iy, y > 0, and 



G(z) = Go(z) + Go(z)T(z)Go(z) 



(13) 



(14) 



with Gq{z) 



z-Ho 



Being equivalent to representation 



Eqs. p3|l and (|14|l express the principle of the su- 
perposition of probability amplitudes. In the ordinary 
quantum mechanics the similar equation plays an im- 
portant role and establishes the connection between the 
evolution operator and the Green operator G{z) which is 
defined by 



G(z) 



1 



z-H' 



(15) 



with H being the total Hamiltonian. Such a form of 
the Green operator follows from the fact that in the 
Hamiltonian formalism the evolution operator satisfies 
the Schrodinger equation. In the canonical formalism 
the T-matrix is defined by Eq. (|14|l starting with the 
Green operator of the form (|15|l . In the formalism of the 
GQD the T-matrix plays a more fundamental role. It is 
defined by Eq. Q, and in this case the starting point is 
representation (Q) with S{t2,ti) being the contribution 
to the evolution operator from the processes in which 
the interaction begins at time ti and ends at time t2. In 
this case the operator G{z) itself is defined by Eq. 114|) 
via the T-matrix. This is a more general definition of the 
Green operator, since representation is a consequence 
of the principle of the superposition of probability ampli- 
tudes and must be valid in any case while the evolution 
operator can be represented in the form (|13|l with the 
operator G{z) given by Eq. (|15|) only in the case where 
the interaction in a system is instantaneous. 

As has been shown in Ref. , there is a one-to-one cor- 
respondence between the character of the dynamics and 
the large-momentum behavior of the T-matrix. If this 
behavior satisfies the requirements of ordinary quantum 
mechanics, then the interaction in a system is instan- 
taneous and the dynamics is Hamiltonian. In the case 
where this behavior is "bad", i.e., does not meet the re- 
quirements of the ordinary quantum mechanics, the in- 
teraction generating the dynamics of the system must 
necessarily be nonlocal-in-time. Let us now illustrate this 
point by using the model developed in Refs. 0, |^ as a 
test model demonstrating the possibility of the extension 
of quantum dynamics provided by the GQD. This model 
describes the evolution of the system of two identical non- 
relativistic particles with the mass m whose interaction 
is separable, and hence the interaction operator has the 
form 

<P2\Htl{r)\p, >=<^*(p2)^(pi)/(r). 



where p is the relative momentum of the particles, and 
/(t) is some function of the duration time r of the in- 
teraction in the system. It is assumed that in the limit 
IpI oo the form factor </?(p) behave as 

'/'(p)-ipr", (iPi->«^). (16) 

In this case the general solution of Eq. JSJ is 

'y'*(P2)V'(Pl) 



(p2|T(z)|pi) = 



ga +[Z-a)j (z-E,)(a-E,) 



.(17) 



where ga = t{a), and a G (— oo,0]. In the case a > i, 
the amplitude (p2|T(z)|pi) given by Eq. (|17|l . tends to 
A(y5*(p2)',5(pi), where 



A - 9a' + 



d^k |y(k)p 
(27r)3 a~Ek 

From Eqs. © and 1)10(1 it follows that in this case the 
interaction operator hI^^I (t) should be of the form 

< P2 1 i/(:i(r) I pi -2zA<5(r)^*(p2)^(pi), 

and hence the dynamics of the system is Hamiltonian 
and is governed by the Schrodinger equation with the 
potential (p2|^|Pi) = ^'P*{P2)'p{Pi)- In the case a < \, 
the T-matrix H17() tends to zero as — >oo 

(p2|T(z)|pi> ^ ^*(P2)^(P1) 

1^1 ^ OO 

X (6i(-z)"-U52(-z)2"-i) +o(|zp"-i), 

where bi — — 47rcos(a7r)m"~i and 62 = 6i|a|^~" — 
62(M(a) + g-i) with 



M{a) = 



-2a 



(27r)3 



a - Ek 



In this case the dynamics is non-Hamiltonian and, as 
it follows from Eqs. lO and H1U|) . is generated by the 
nonlocal-in-time interaction operator 



< P2 1 i/i:)(r) I pi > = ^*(P2)^(P1) 



X aiT " 2 + 



-2a 



(18) 



where ai = Airi cos{a7T)m"~^T~^{^—a) exp[z(— ■|-|-l;)7r], 
and 02 is a free parameter of the theory. The solution of 
Eq. © with the interaction operator (|18|l is of the form 



(P2|T(z)|pi) 



fef</7*(p2)(^(pi) 



-&2 + &i(-z)^-" -M{z)hl 



III. THE PIONLESS THEORY OF NUCLEAR 
FORCES 

The Weinberg program for low energy nucleon physics 
employs the analysis of time-ordered diagrams for the 2N 
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T-matrix in ChPT to derive a NN potential and then to 
use it in the LS equation for constructing the full NN 
T-matrix. Obviously the starting point for this program 
is the assumption that in the nonrelativistic limit ChPT 
leads to low energy nucleon dynamics which is Hamilto- 
nian and is governed by the Schrodinger equation. How- 
ever, the fact that the chiral potentials constructed in 
this way are singular and lead to UV divergences means 
that this is not the case. At the same time, the GQD 
allows one to analyze the predictions of ChPT without 
making a priori assumptions about the character of low 
energy nucleon dynamics: This character should results 
from the analysis. Let us consider, for example, the low 
energy predictions of ChPT for the 2N system in the 
^5*0 channel. At very low energy even the pion field can 
be integrated out, and the diagrams of the ChPT take 
the form of the diagrams being produced by the effective 
Lagrangian containing only contact interactions among 
nucleons and derivatives thereof JJJ. From the analysis 
of time-ordered diagrams of this theory it follows that 
the 2N T-matrix in the ^S'o channel must be of the form 

oo 

(p2|T(z)|pi) = J2 PTpl'^tnmiz), Ipil < A, IP2I < A, 

n,m— 

(19) 

where p^ is relative momentum of nucleons, and the 
functions t„™(z) are of order 0{(Q/A)2"+2™}. Qn the 
other hand, the 2N T-matrix must satisfy the GDE 
of the form ©. From Eq. (fT^ it follows that for 
T(z, p2,pi) = (p2|T(z)|pi), we can write 



T(z,p2,pi) = ^p2«^,(z,pi), IP2I < A, 

OG 

T(z,p2,Pi) = ^p?"t„(z,p2), Ipil < A, 

00 

Tn{z,p)=Y,P^""tnmiz), |p|<A, 
00 

fn{z, P) = ^ P^"'UUz), IpI < A. (20) 



m=0 

By using these expansions, Eq. ((SJ for the 2N T-matrix 
can be written in the form 

o?T(z,p2,Pi) ^, , , sdMoiz) 
= ro(z,p2)To(z,pi) — 

az az 

+ ^(z,P2,Pi), (21) 



with 



. (Pk 1 m 

M^){z) = zm I 7777:7^7 , ^ = —V—zm, 



J'{z,P2,Pl) 



(27r)3 {zm - k^)Ek At: 



(Pk Ti(z,p2,Pi,k/A) 
(27r)3 (z - 



(22) 



.Fi(z,p2,Pi,k/A) = T(z,p2,k)T(z,k,pi) 
- To(z,p2)To(z,pi). 

The function .F(z, p2,Pi) is of order 0{Q/A}. This is 
because, Eq. (|22(l can be rewritten in the form 



•^(2:,P2,Pl) 



m f (fik Ti(z, p2, Pi, k) 



A 



(27r)3 (r7iz/A2-fc2)2' 



where k = k/A and Ti (z, p2, Pi, k) = ojfc^}^ k ^ 0, due 
to the above substraction. This allows one to construct 
the 2A^ T-matrix in the spirit of the effective theory by 
expanding it in powers of Q/A 

T(z, p2, Pi) - t(2A^)(z, P2, pi) + O {(Q/A)2-^+i} , 

where T^^-^) (z, p2, pi) is the sum of the first J\f terms in 
the expansion. The function T^'^^\z, p2,Pi) must sat- 
isfy Eq. lEU with an accuracy of order O {(Q/ A) 2^+1} 



dT(2^)(z,p2,Pl) 



dz 



dMojz) (2^) 



7^(2^) r, ^'r(2A/')/ n 
-'0 {Z,P2)-'q [Z,pi) 



^(^^)(z,p2,Pi) {1 + 0{{Q/Af^+'}) 



dz 
where 



^'-'(.,P.P.)^-/^^(r<-i(.,P.« 
xT<»«(.,k,pi) - r;»''(.,p,)r;»'>(.,p,) 



Taking into account that T(z,p2,Pi) is of order 
0{Q/A}, Eq. (|23(l can be rewritten in the form 



dT(2A^)(z,P2,Pl) 



dz 



dMojz) 
dz 



fr\z,P.)rr\z,p,) 



+ .F(2(A^-1»(Z,P2,P1)" 

X (l + 0{(Q/A)2-^+i}) . (23) 



This equation allows one to construct the 27V T-matrix 
order by order. 

At leading order Eq. (|23|) takes the form 

'"^'"'^^f'P^^ = rr(.,P.)rr(.,Pi) 



dMa{z) 
dz 



(l + 0{g/A}). (24) 



By using this equation and the above expansion of 
T(z, p2,pi), for tQo{z) at leading order (LO) we get the 
equation 



dC{z) 
dz 



{Ci4'-^{l + 0{Q/A} 



The solution of this equation is 

4o (^) - (C(7' - Mo{z))-' (1 + 0{g/A}), (25) 
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where Co is some constant. Correspondingly, for 
T^^\z, p) we have the equation 



rf7o^"^(^.P) _ ,(0) 



dz 

whose solution is 



4?(^)ro(^^,p)^(i + o{g/A} 



(0) 



(z, p) = 4°^ (z)*o(p) (l + 0{Q/K\) , (26) 



where ^o(p) is some function that for |p| < A can be 
expanded as 



*o(p) = E^2„/", |P|<A, 



(27) 



11=0 



with Co — 1. In the same way, for T^*-"^ (z, p), we get 

f^'\z,Y>) = 4o (z)*o(p) (1 + 0{Q/A}) , (28) 
where the function ^o(p) satisfies the condition 

oo 

^'o(p) = E4y", |P|<A. 

Tl = 

Now, substituting 7^'"^ (z, p) and TJ,'"^ (z, p) into Eq. (Pl|l . 
we can obtain the amplitude r(z, p2, Pi) at LO 



r(o)(- 



^^,P2,PlJ 



*o(P2)*o(Pi)(C^o"'-A^o(^)) 
*o(P2,Pi)l (1 + 0{Q/A}), (29) 



where $o(P2,Pi) is an arbitrary function satisfying the 
following conditions 

oo 

$0(Pl,P2) = 5Zp?"'^2„(P2), |Pi|<A, 
n=l 
oo 

*o(P2,Pl) = 5I?'?"^2n(P2), |Pi|<A, 

n=l 
oo 

<^2n(p) = X! C2n,2mP^", |p| < A, 
m— 1 

oo 

<^2«(P) = ^ C2™,2„P^'", IpI < A. 



m— 1 



These conditions result from the fact that, for the func- 
tion $o(P2, Pi) to be consistent with the structure of the 
2N T-matrix shown in Eq. p9|l , it should be of the form 



2m 
Pi ) 



*o(P2,Pl)= XI C2„,2mP2"?: 
n,m— 

IpiI < A, |p2| < A. 



(30) 



On the other hand, Co, 2m and C2n,o should be zero be- 
cause if this is not the case the terms Co,o and Co.2nPi" 



have to be included into ioo^(z) and Tq^^(z,pi), respec- 
tively. As has been shown in Ref . 3 , for the solution of 
Eq. to be unitary, it must satisfy the condition 

(p2|T+(z)|pi> = (p2|T(z)|pi), z G (-00,0). 
From this and Eq. it follows that 3'o(p) = *o(p)' 

<^2n(p) = ¥'2ri(P)> C2„,2m = C2„i^2n: and $o(P2,Pl) = 
^'o(Pl,P2). 

The fact that Eq. (j^HJl represents the 2A^ T-matrix 
at leading order means that at extreme low energies the 
true 2N T-matrix in the ^5*0 channel is described by this 
equation. In other words, at z = the true 2A^ T-matrix 
should coincide with T('^)(0, P2, Pi) given by Eq. iP^. 
and hence the higher order corrections must not change 
(p2|T(z)|pi) at z = 0. This puts the boundary condition 
onT(2A^)(z,p2,Pi) 

T(2^)(0,P2,Pl) = *S(P2)*o(Pl)Co + $o(P2,Pl). (31) 

By using Eq. H23|l with this boundary condition, one 
can obtain the 2A^ T-matrix at any order. However, the 
functions ^'o(p) and $o(P2!Pi) that occur the bound- 
ary condition (|31|l arc not known exactly because such 
a knowledge implies that we know the details of the un- 
derlying dynamics (the relevant domain of the definition 
of the functions 4'o(p) and <i>o(p2,Pi) spreads to high 
energies). As we show below, for obtaining T(z, p2,Pi) 
with an accuracy of order 0{{Q/KY^+^} it is sufficient 
to know a few constants C2n and Cn,m that appear in the 
expansions (|27ll and H30|l . and a few parameters charac- 
terizing integral properties of the functions 5'o(p) and 
<i>o(p2,Pi)- This allows one to describe low energy nu- 
cleon dynamics in terms of a few low energy constants 
(LEC's) without any knowledge of the details of the dy- 
namics at high energies. 

For describing low energy nucleon dynamics at leading 
order it is sufficient to restrict oneself to the first term 
in the expansion (|27|l of the function \l/o(p), while the 
function $o(P2,Pi) can be omitted because Co,o = 0. In 
this case the LO 2A^ T-matrix takes the form 



(p2|T(°)(z)|pi) = (Co-i-£V^ 



(32) 



where the low energy constant (LEG) Cq parametrizes 
the effects of high energy physics on low energy nucleon 
dynamics. This constant can be obtained by fitting to 
the 2A^ scattering data. Since the scattering amplitude 
Mp) = -(p|T(°H^^p+«0)|p) at p = should be described 
by the LO 2A^ T-matrix, the constant Co is related to the 
scattering length a as Co = A-Ka/m. Note that the stan- 
dard EFT approach yields the same expression for the 
LO T-matrix 0, [H |li|. It is remarkable that the T- 
matrix H32(l . which in the EFT approach is obtained by 
performing regularization and renormalization of the so- 
lution of the Schrodinger (LS) equation, is not a solution 
of this equation. At the same timej_it is a solution of Eq. 
^ with the interaction operator 



(P2 



5(")(z)|pi) 



47r 



167r2 



m^—zm CQirfizm 



(33) 
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The corresponding generalized interaction operator 
^int {''') is of f tie form 



P2i4:)Mipi> 



27r\3/2i 
m 



TT V 



/2 



(34) 



where 9 = (7^)'^^^C'o'^(l + This operator is nonlocal 
in time and this is the only reason why in this case the 
GDE cannot be reduced to the LS equation, and hence 
the T-matrix is not its solution. 

At the same time, for using r'^°)(z, P2, Pi) in Eq. H23|) 

(more precisely the first term Tq'^\z,p) in its momen- 
tum expansion l|2()(l '). one cannot restrict oneself to the T- 
matrix shown in Eq. H32|l . This is because in this case the 
LO T-matrix is used for obtaining higher order correc- 
tions to T(z, p2, pi), and hence must contain a more de- 
tailed information than the T-matrix responsible for de- 
scribing nucleon dynamics at leading order. Thus, while 
by definition T(z, p2,pi) is equivalent to (p2|r(2;)|pi), 
the amplitudes T^^^\z, p2, pi) that are used in Eq. 
and (p2|r''^^''(z)|pi) that control the nucleon dynamics 
with an accuracy of order 0{{Q/A)^^+^} represent dif- 
ferent approaches to the 2N T-matrix: In contrast with 
(P2|T(2-^)(z)|pi) the amplitude T'^^^^z, p2,Pi) must 



contain the information that allows one not only to de- 
scribe low energy nucleon dynamics at leading order but 
also to obtain the higher order corrections. At LO all the 
needed information is contained in the functions "00 (p) 
and <i>o(P2,Pi) that appear in Eq. (|29() . In principle 
they could be derived from QCD. However, for describ- 
ing low energy nucleon dynamics at LO one need not to 
know these functions exactly. It is sufficient to know the 
first terms in their low energy expansions. This manifests 
itself in the fact that the T-matrix describing low energy 
nucleon dynamics at LO is momentum independent. At 
the same time, as we show below, for obtaining higher 
order corrections a better knowledge of the properties of 
the functions 4>q{p) and $o(P2,Pi) is needed. 

At next to leading order {NLO) Eq. if^ reads 



dT(2)(z,P2,Pi) 



dz 



rJ^)(z,P2)ro^^'(z,Pi) 



(2), 



x^^+^(")(z,P2,Pi)] (1 + 0{QVA^}), (35) 



and hence ^qo (-^) satisfies the equation 



dz 



dz 



(1 + 0{QVA^}), 



where M(z) = zm J (^Ji^l%%^ with Ti(fc/A) = |*o(k)P - 1. The fmiction M{z) can be expanded as 



^ / Ti(fc/A) , ,2 
M{z) = zm / - — ^ / — = (zmY 



d^k Ti(fc/A) 



d^k Ti(fc/A) 



(27r)3 {zm ~ k^)Ek ' ' J (27r)3 {zm - k^)Ekk^ J {2tt)^ k^ 

oo / n \ oo 

= X! i^'^T Mo{z)^C*2,C2(ri-z) - J2n{zmY = ^ M2n{z), 



i=0 



n=l 



(36) 



(37) 



i=0 



with 

n 

M2n{z) (zm)" Mo{z) ^ C*2^C2{n-i) ~ J2n 

where 

= mA 



d^k F„{k/A) 
(27r)3 fc2(»+i) 

-2n+l f <P~k F^k) 



(27r)3 k2(n+l) 



(38) 



(39) 



F„(fc/A) = |vI/o(k)|2 - ^ ^C^,C2(,_,)fc2^ 

i=0 j=0 

and k = k/A. From this it follows that the function 
M2„(z) is of order OjQ^^/A^"}, and hence in Eq. ^ 



it is sufficient to keep only the first term in the expansion 
(|37|l . Thus Eq. (|36|l can be rewritten in the form 



dz 



dMo{z) / (0) 



dz 



dM2{z) 



dz 



X (i + o{qVa^}) 



(40) 



The boundary condition for this equation should be as 
follows 



4^(0) = Co. 



(41) 



This is because at z = the exact 2A^ T-matrix should 
coincide with the LO T-matrix and this puts the bound- 
ary condition on T(^^^(z, P2, Pi). The solution of Eq. 
(|40|l with this boundary condition is 



tg\z) ^ [C^' - Aloiz) - M2iz)\ 
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Now we can obtain the amplitude Tq{z,p) at NLO that, 
as it foUows from Eqs. it^ . (j^ . and should satisfy 
the equation 



dz 



where 



A4{z, p) = zm 



X (1 + 0{qVA^}) , (42) 



(27r)3 {zm-k^)Ek ' 
The function A^(z,p) can be expanded as 

M{z,p) = M2{z,p) {l + O {QVA'}) 

AT 

-E-^2„(^,p)(l + 0{(Q/A)2-^+^}) 



with 



■A^2n(2:,P) = (zm)"Mo(z)^C2j(/72(«-o(P) 



j=0 



(zm)"a2„(p) 



where 

a2n(p) = ?™ 



(27r)3 



„-2(n+l) 



^'o(k)$o(k,p) 



n-l i-1 

-H^^^H'^2i¥'2(i-»)(p) 

j=l i=0 

This allows one to rewrite Eq. (|42l) in the form 



d%^'\z,p) 



dz 



(4°'(z))%o(p) 



4^o^(-)ro^^^(.,p)^ 

dM2(z) , ^(0)^ ^d.M2(z,p) 



dz 



dz 



X (1 + 0{qVA^}), (43) 



From Eq. (EU it follows that at z = 

%^'''\0,p) = Cmoip) = CoMp)- (44) 
Solving Eq. H43|) with this boundary condition yields 
„^ _ «'o(p) + Al2(2:,p) 



V'{z,p) 



C'o^ ~ Moiz) - Ahiz) 

X (i + o{qVa'}). 



(45) 



" (2) 

In the same way, for 7^^ ^(z,p), we get 



tj'^(^,p) 



V|/^(P)+.M2(Z,P) 

Cq-i -Mo(z)-M2(z) 
x(l + 0{OVA^}), 



(46) 



where 7^2(2, p) = Ml{z* ,p). Substituting these expres- 
sions into Eq. (|35|l . we get the equation 



dT(^)(z,P2,pi) 
dz 



*S(P2)*0(P1)(4?(^))' 



(*S(P2)-M2(Z,P1)+^2(Z,P2)*0(P1)) {tg\z) 



fdMoiz) dM2{z) 
\ dz dz 

(2), ,\^dMoiz) 



dz 



,+(0)/ \ Lr,*t .dM2iz,Pi) , dM2{z,P2) 
+*oo (2) *o(P2) — + *o(Pi)- 



dz dz 
Solving this equation with the boundary condition H31|l yields 

■*?i(P2)^'o(Pl) + *5(P2)-A^2(^, Pi) + A^2(2, P2)^'o(Pl) 



i + o{qVa^} 



(47) 



r(2)(z,p2,pi)- 



Mo{z)^M2iz) 



$o(p2,Pi) (i + 0{QVA'}). (48) 



r 



In describing low energy nucleon dynamics at NLO the 
terms containing the functions ^'o(P2jPi)j f2{p) and 
a2(p) can be omitted, and hence, for (p2|T(^^(2:)|pi), 
we get 



(P2 



\T'^'\z)\Pi) 



^5(P2)'I'o(Pl) 

~ Mo{z) - M2{z) 

X (i + o{qVa'}). 



(49) 



For momenta below A it is sufficient to keep the first two 
terms in the expansion (|27|l . and the NLO T-matrix can 
be represented in the form 

{P2\T^'H^)\PI) 

{1 + cIpI){1 + C2pI} 



Cq^ - Moiz)(^l + 2Rec2zrnj + J2zm 
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X (l + 0{gVA3}), |pi|<A, |p2|<A. (50) 

Thus, at NLO the effects of the high energy physics that 
in the full theory are described by the function 5*0 (p) are 
parametrized by the two LEC's C2 and J2- The constant 
C2 determines the second term in the expansion H27|l of 
the function 4'o(p), while characterizes its integral 
properties shown in Eq. (|39|l . Really for given C2 and 
J2 Eq. H49|l represents the set (we will denote this set 
^2) of the solutions of the GDE that coincides with the 
true 27V T-matrix with an accuracy of order OjQ^/A'^} 
provided that C2 and J2 are the LEC's chosen by nature. 
Equation (|29|l also represents a solution set of the GDE. 
The LO T-matrix given by Eq. I|32|l belongs to this set, 
and hence can be used for describing nucleon dynamics 



at leading order. In contrast, Eq. H50|l does not represent 
the T-matrix belonging to the set ^2 because it does not 
represent the solution of the GDE at all, and can describe 
the LO 2A^ T-matrix only for momenta below the scale A. 
For this reason, to describe nucleon dynamics one should 
use Eq. (|49|l with the function ^'o(p) satisfying the two 
conditions: it must have the low momentum expansion 
(I27|l where second term is determined by the LEG C2, 
and is related to the LEG J2 by Eq. These LEG's 

should be derived from low energy experiment, and first 
of all from the NN scattering data. 

The scattering amplitude that results from the NLO 
T-matrix shown in Eq. H50|) is 



= _(p|T(2)(z = i?p + zO)|p) = - 
1 



{\ + e2v'){\^c2v') 



An 
m 



Cq^ - Mo(z)(l + 2Rec2zm) + J2zm 

+ o{qVa'} 



(i + o{qVa^} 



47r2C^'Rec2- J2 2 
■ p 



ip}^ 



(51) 



Taking into account that for p << A this equation must 
reproduce the effective range expansion 

A(p) = _ zap + (aro/2 - a')p^ + 0{p'}], 

m 

we get the relations between our LEG's and scattering 
length a and the effective range tq 



Co 



47ra 
m ' 



2Rec, 



47ra 



-J2 



2 



(52) 



Thus, at NLO we have two equations that relate the 
three LEG's Co, C2 and J2 to the scattering data. This 
means that the two nucleon scattering data are not suffi- 
cient to fix these LEG's, and fitting to experimental data 
where the off-shell properties of the 27V T-matrix mani- 
fest themselves is needed. For example, one can use the 
pp bremsstrahlung observables that are very sensitive on 
the low energy NN interaction. 



As we noted, for constructing the T-matrix at higher 
orders we have to deal with the amplitude T^'^\z, P2, Pi) 
given by Eq. 148|) . This amplitude contains the functions 
*J'o(P2,Pi), V2(p) and ^'2(p) that do not manifest them- 
selves at NLO. These functions should come into play at 
higher orders. At next-to-next-to-leading order (NNLO) 
Eq. takes the form 



dT(4)(z,p2,pi) 



dz 



T^'\z,p2)Tr{z.l>i) 



(4)/ 



dMo(z) 



+^(2)(z,P2,Pl) 



dz 

o{qVa'}). 



Gorrespondingly, by using Eqs. 145|l . (|4(i(l . and 149|l for 



tpQ^ (z) we get the equation 



4o(^) _ ,\^ dMo{z) , ^^(2), j dM2{z) , dMijz) 



d^k ^*oik)M2iz,k) + ^oi^)M2{z,k) 



X (i + o{qVa^}) = [(4o'W 



(2^)3 (z - 

+ (4o W)' ^{^^2(^) + M4(z) -f C2,2(^m)2M2(z) 



dMp{z) 



where 



AA — m 



d^k ^o(k)^;(k) 
(27r)3 fc4 



BA—m 



d^k ^g(k)Q2(k) 

(27r)3 k^ 



(53) 
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Solving this equation with the boundary condition (|41() yields 

,(4). . (i + o{gVAn) 



1 + 2Rec2Zin 



leap + 2Re(c4 - A4) + C2,2Mo(z) (zm)2| + J^zm + 



Ji - ReB4 



This expression for t\^Q{z) can than be used in the equation for T^'^\z,'p) 



d^k *S(k)X2(^,k)- 



d^k *o(k)A^2(^,k) (2), , f d^k [*o(k)+A^2(^,k)]$o(k,p 

'-00 



(27r)3 {z-Ekf ' (27r)3 

The solution of this equation with the boundary condition 144f) is 



(z - Ek Y 



(2^)3 (z - Ek^ 

i + o{qVa^} 



(4) 



where 



(z,p) = 4'o'W(*o(p) + Al2(^,p) + Al4(^,p) - ^(z™)2Mo(z)04(p) + l{zmfl4{p)) (l + 0{QVA'}) , 



d^k <^2(k)$o(k,p) 



fc4 



, 74 (p) = m 



(27r)3 A:4 



In the same way for 7'q^\z,p) we get 

"(4)^^ „\ _ .(4) 



7o^'^(^,p) = 4o^(z)(*S(p) +^2(^,P) +^4(^,P) - ^(zm)2Mo(z)0:(p) + i(zm)^4*(p)) (l + O {QVA^}) . 



Now we can rewrite the equation for T^*^' (2:, p2, pi) in the form 



dr(4)(z,p2,pi) 



dz 



rr(.,P2)rr(.,Pr)^--^^^ 



d^k *5(k)A^2(^,k) + *o(k)^2(2,k)^ ^(2), , 



ro^'^(^,Pi) 



^3;. $o(P2,k)(«'5(k) +7^2(2, k) 



dM2(2) dAf4(z) 
^3;. (^-olk) +X2(2,k))$o(k,Pi) 



r„(^)(z,P2)To(^)(z,Pi)( 



(27r)3 {z - Ekf dz 

Solving this equation with the boundary condition (|31f) yields 



(2^)3 (z - EkY 

d^k $o(P2,k)$o(k,Pi) 



(2^)3 (z - Ek)k-' 



(i + o{gVA5} 



TW(z, P2,pi) = rJ")(z,p2)Tj"^(z,pi) tf^iz) +$o(P2,Pi)-2™T4(P2,Pi 



(54) 



where T4(p2,pi) — m J i^^J)^ . In describing low energy nucleon dynamics the terms including the 

function A^2(^jP) can be neglected, and hence for (p2|T('*^(z)|pi) we can write 



(P2 



|rW(z)| 



pi; 



{*S(P2) + >l2(z, P2)}{*o(Pi) + -M2(z, Pi)}4o^(z) + a>o(p2, Pi)] (1 + OiQ^A^}). (55) 



For momenta below A this equation can be rewritten in the form 



(P2|r('*)(z)|pi) = (l + c*2pI + clpj + Cl2zmMo{z)pj - Alzmpfj (l + C2pl + Cipf + C2,2zmMoiz)pl - Aizmpf^ 



X 4Hz) + C2,2P^p?](l + 0{QVAn) 



(56) 



Here we have used the fact that for |p| < A the function 02 (p) can be expanded as 

a2{p) = Aap^ + A&p^ + |p|<A. 
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By using this equation and Eq. for the NNLO 2N 

scattering amplitude we can write 



in 



m a ^ 



„4 



-p 



a-' H 
ri/2 



„4 



where n = • X with X = Cp"^ 



|C2| 



- 2Re(c4 - Aa) - 



(2Rec2)2 + Co-T2,2j - + ReS4 + 2J2Rec2. This is 

the Kaplan-Savage- Wise (KSW) expansion of the effec- 
tive theory [T^. At the same time, Eq. (|56|l allows one 
to obtain the more general expansion 



A{p) = , . 

{C2{^Ji)p^fnl{^l + ip)/i^: 



C2W 



Ci{li)p'' 



Co(jj.)r. 



4ir 



■ip + ipW 



with CqI^j,) 



X 



Co{fJ.)r 
iir 



C2(m) 



,(57) 



:72jC2(/x) and C4(m) - - • -0 
where fj, is an arbitrary parameter. 

We have demonstrated that, by using Eq. (|^ and 
taking into account that the 2N T- matrix in the ^Sq 
channel should be of the form H19() . one can construct 
this T-matrix order by order in the expansion in pow- 
ers of Q/A. For example, the 2N T-matrix up to NLO 
and NNLO is given by Eqs. and respectively. 
It should be noted that the way in which the above ex- 
pansion of the 2N T-matrix has been obtained is not a 
regularization procedure. Moreover, the key point of our 
approach is that the 2N T-matrix must necessarily sat- 
isfy the GDE without regularization. This requirement 
directly follows from the fact that at low energies nu- 
cleons emerge as the only effective degrees of freedom, 
i.e., the probability to find a system, the initially was in 
a two-nucleon state, in a state including other degrees 
of freedom is negligible. Indeed, the above means that 
the evolution operator defined on the 27V subspace of the 
Hilbert space should be unitary. From this it follows im- 
mediately that the 2N T-matrix must satisfy the GDE. 
This is because this equation is an inevitable consequence 
of the requirement that the evolution operator is of the 
form and satisfies the unitary condition. Our expan- 
sion of the 2N T-matrix is based on the fact that the 
GDE must be satisfied at any order of the expansion in 
powers of Q/A. Equation (|23|) represents the approxi- 
mation to the GDE that allows one to obtain the 2N 
T-matrix with an accuracy of order 0{{Q/A)^^+^}, if it 
is known with an accuracy of order 0{{Q/Ay^-^}. As 
we have shown, in this way we get the same expansion 
of the scattering amplitude that in the standard pionless 
EFT is derived by summing the bubble diagrams and 
performing regularization and renormalization. The ad- 
vantage of the above way of the expansion of the effective 



theory of nuclear forces is that it allows one to construct 
not only the scattering amplitude but also the off-shell 
T-matrix. As we show in the next section, this provides 
a new way to formulate the effective theory of nuclear 
forces. 



IV. THE EFFECTIVE THEORY OF NUCLEAR 
FORCES. 

As follows from Eq. in order to describe the low 
energy dynamics we have to obtain the T-matrix for any 
z relevant for the low energy theory. Equation ^ allows 
one to obtain the operator T(z) for any z provided that 
some boundary condition on this equation is specified. 
The boundary condition © means that the most of the 
contribution to the operator T(z) in the limit \z\ — > 00 
comes from the processes associated with the fundamen- 
tal interaction in the system under study that are de- 
scribed by the operator B{z). On the other hand, for the 
low energy theory to be consistent, the operator B{z) 
must be determined in terms of low energy degrees of 
freedom. This means that, despite the boundary condi- 
tion formally implies that z must be let to infinity, 
really one has to restrict oneself to a " high" energy region 
which is much above the low energy scale but is well be- 
low the scale of the underlying high energy physics. We 
will use V to denote this "high" energy region. The fact 
that the energy region V lies much above the scale of the 
low energy dynamics implies that at such "infinite" en- 
ergies the main part of the contribution to the operator 
T{z) comes from processes that can be thought of as a 
"fundamental" interaction in the low energy theory, and, 
as a result, the interaction operator B{z) is close enough 
to the true 2N T-matrix. Thus, in order to take into 
account the fact that any theory has a range of validity, 
instead of the condition , we have to use the following 
boundary condition on Eq. (|SJ): 

(^-21^2)1^1) = {i:2\B{z)\^/j^)+0{h{z)},zeV{58) 

In order for Eq. JHJ with the boundary condition f5H|) to 
have a unique solution, the interaction operator must be 
close enough inside domain V to the true 2N T-matrix, 
i.e., the function h(z) must be sufficiently small for z g 
v. This means that the operator B{z) must satisfy the 
equation 



d(V^2|g(^)|^l 

dz 



E 



{^2\B{z)\n){n\B{z)\iji 



[z - i?„)2 
+ 0{\z\-^h{z)] , zeV. (59) 



In the low energy theory of nuclear forces all processes 
that in ChPT are described by irreducible diagrams in- 
volving only two external nucleons can be considered as 
such " fundamental" processes. Here irreducible diagrams 
are 2A'^ irreducible: Any intermediate state contains at 
least one pion or isobar. It is natural to expect that the 
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relative contribution of reducible 2N diagrams tends to 
zero as z increases, and, in a region that lies much above 
the scale of low energy nucleon dynamics, the main part 
of the contribution comes from the "fundamental" pro- 
cesses that are described by the 2N irreducible diagrams, 
and therefore the interaction operator B[z) should de- 
scribe the contributions from these processes. On the 
other hand, because of the separation of scales provided 
by QCD the above "high" energy region of the low en- 
ergy theory lies still much below the scale of the underly- 
ing physics. In other words, the amplitudes {ip2\B{z)\il]i) 
that constitute the interaction operator and generate low 
energy nucleon dynamics really are low energy (in the 
scale of the underlying theory) amplitudes that in ChPT 
are described by the irreducible 2N diagrams. In prin- 
ciple they can be obtained within the underlying high 
energy theory in some low energy limit, and can be di- 
rectly used as building blocks for constructing the low en- 
ergy theory. Thus the GQD allows one to build a bridge 
between QCD and low energy nucleon dynamics. It is 
hoped that in the future it will be possible to obtain 
the amplitudes {'4)2\B{z)\il}i) in terms of QCD with such 



accuracy that the corresponding operator B{z) will de- 
termine a unique low energy 2N T-matrix. There is no 
reason to believe that this operator must necessarily gen- 
erate the Hamiltonian low energy dynamics. In fact, as 
it follows from Eq. (|11|1 . for this (p2|-B(2:)|pi) must have 
a negligible dependence on z inside the domain V. How- 
ever, as we show below, such a behavior of (p2|-B(z)|pi) 
is at variance with the symmetries of QCD. On the other 
hand, in the region V QCD must reproduce the 2N T- 
matrix that satisfies the GDE and hence can be used as a 
NN interaction operator governing low energy dynamics. 

In the case of the pionless theory the problem of con- 
structing the generalized interaction operator actually is 
a special case of the problem of obtaining the 2N T- 
matrix when we restrict ourselves to the domain T). We 
should also keep in mind that, as it follows from Eq. 
(|^. in the domain V the operator B{z) coincides with 
the true T-matrix with an accuracy of order 0{h{z)}. 
In other words, the operator B),g'j{z) that describes the 
NN interaction up to next-to-leading order is related to 
the operator T'^'^\z) in the domain V by 



P2|T(^^W|pi> 



*5(P2)*0(P1) 



l + 0{QVAn =-*S(P2)*o(Pi) 



Co~Mn{z)^M2{z) 

M^\z){l - 2Rec2ZTO) + M^''{z){c^^ + J^z-n^ (l + 0{M^^ {z)C^'^}) (l + 0{QVA^}) 

= (P2|i?g^/2)|pi) (1 + 0{Afo-3(z)Co-2}) , z e 



r 



with 



(P2|i?ffU^)|Pl> 



'e//V-;iHi/ - -*S(P2)1'o(Pl) 

Mo"i(z)(l - 2Rec2zm) + M^'^{z)(g^^ + J^azm^ 

X (1 + 0{gVA3})(fiO) 

In this equation the effects of the high energy physics are 
parametrized by the EEC's Co, C2, and J^- These EEC's 
put the constraints on the form factor ^o(p): The LEC 
C2 determines its low momentum properties of the form 
factor (p) with an accuracy of order (Q/A)'^, while the 
constant J2 determines integral properties of this func- 
tion relevant for calculations at the same order. The 
above means that, in contrast with the operator B(£) 
which unequally determines the relevant solution of Eq. 
0, the effective interaction operator B'f^^iz) determines 
the set 572 of the solutions of this equation that coincide 
with the true 2N T-matrix with an accuracy of order 
0{((3/A)^}. From the above expansion of T(^)(z, P2, Pi) 
it follows that the function h(z) which determines the 
size of the terms in the interaction operator that can be 
omitted should be of order 



h{z) 



M^\z)\}. 



(61) 



The effective interaction operator B\^j{z) given by Eq. 
(|60|l actually represents a set of the interaction operators 
that generate the T-matrices belonging to the solution set 
5I2 • Each of these operators corresponds to some definite 
function ^'o(p) satisfying the above requirement. Such 
an operator is of the form 



(p2|B(z)|pi 



-M'S(p2)*o(Pi)[A/o"X^) 



-M^\z)C^ 



and generates the T-matrix 



(p2|T(z)|pi) ^ 



«'S(P2)*0(Pl) 



C^' - Mo{z) - M{z) 



(62) 



(63) 



where M{z) is given by Eq. (|37|l . It should be noted 
that in contrast with the effective interaction B),^^{z) 

the operator H62|l contains only two terms Mq^{z) and 
Mq'^{z)Cq'^ that determine its dependence on z. This is 
because the form factor vE'o(p) that occurs in Eq. ll^ 
is assumed to be known, and the operator H62|l is close 
enough in the domain T) to some separable T-matrix with 
this form factor to unequally determine it. The function 
M{z) describing the dependence of this T-matrix on z is 
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completely determined by the form factor (p). From 
the fact that the form factor satisfies the above require- 
ments and hence the T-matrix f^K^ belongs to the so- 
lution set ^2 corresponding to given LEC's C2 and 
it follows that the first term in the expansion of M{z) 
shown in Eq. IpTfjl 

M(z) = [2Rec2zmMo{z) - J2zm] (l + 0{QVA^}) 

corresponds to the same constants. In contrast, the ef- 
fective interaction operator B^^j(z) determines the solu- 
tion set U,2 as a whole not the solution corresponding to 
a given form factor, and the dependence of this operator 
on z that is described by the term zmJ^2 — '^^^C2zmMQ{z) 
puts the above constraints on the function ^'o(p). 



Equation (|55|l represents the set Q.^ of the solutions of 
the GDE that coincide with the true 2N T-matrix with 
an accuracy of order OjQ^/A^}. This solution set is de- 
termined by the values of the LEC's Co, C2, C4, C2,2j J^2i 
J4, Ai and B^. These LEC's put the constraints on the 
functions ^'o(p) and $o(P2,Pi) that determine the defi- 
nite solutions belonging to the set f24: The constants C2, 
C4 and C2,2 that occur in their low momentum expansions 
shown in Eqs. H27|l and H30() have the same values as the 
above LEC's, and these functions have the integral prop- 
erties that are parametrized via Eqs. H39|l and H53|) by 
the LEC's J2, Ja, A4 and B4. The effective interaction 
operator that generates the solution set fii of the GDE 
is of the form 



{p2\Bi%{z)\p^) 
1 — 2Rec2zm — 



:{m,-\z)[ 



+Mo^{z)(c^^ +J2zm + 



|C2| 



- {«'S(P2) +Al2(z,P2)}{«'o(Pl) +X2(2,Pl)} 

2Re(c4 - ^^4) - (2Rec2)^ + C2,2Mo(z)j (zm)^ 



Ji - ReB4 



(zm)2) } + $o(p2, pi)] (1 + 0{QVA'} 



(64) 



This equation represents the set of the interaction oper- 
ators that generate the solutions of the GDE belonging 
to the set fl^, and for calculations with an accuracy of 
order 0{Q^ / A^} one can use the operator 1)64(1 with any 
functions "00 (p) and $o(p2,Pi) satisfying the above con- 
ditions with the corresponding LEC's. In the same way 
we can obtain the set il.2^f of the solutions that coincides 
with the true 2N T-matrix with an accuracy of order 
0{((5/A)^^+-^} for any Af, and as Af increases one ap- 
proaches to the true 2N T-matrix closer and closer. 

As we have noted the formalism of the GQD prescribes 
to start with the analysis of the theory at high energies 
where in general the structure of the theory is most sim- 
ple, but in the pionless theory this is not the case. In this 
theory the Weinberg analysis of the diagrams in ChPT 
yields the same expansion (|19|l of the 2N T-matrix for 
all relevant energies, and one can use it for obtaining 
this T-matrix not only in the domain 2? but also at low 
energies. However, in the case when pions must be in- 
cluded as explicit degrees of freedom the structure of the 
2N T-matrix can be derived directly from the analysis of 
the diagrams in ChPT only at "high" energies belonging 
to the domain V, where the theory has the most simple 
structure. Indeed, in this case the T-matrix given by Eq. 
(|58|l should be considered as the contribution from all 
relevant contact diagrams in ChPT, and represents not 
the full 2N T-matrix but only its short-range compo- 
nent. The long-range component should be obtained by 
iterating the pion exchange potential that parametrizes 
the long-range part of the NN interaction. However, 
only at high energies the full 27V T-matrix can be repre- 



sented as a sum of these contributions. This point will 
be demonstrated in the next section. 



V. A NEW LOOK AT THE WEINBERG 
PROGRAM 

The Weinberg proposal was based on the assumption 
that the only equation that can govern low energy nu- 
cleon dynamics is the LS (the Schrodinger) equation, and 
hence what one has to derive from the analysis of dia- 
grams in ChPT is an NN potential. However, there is 
no reason to consider that low energy nucleon dynamics 
is necessarily governed by the Schrodinger equation. In 
principle this dynamics may be governed by the GDE 
with a nonlocal-in-time interaction operator when this 
equation is not equivalent to the Schrodinger equation. 
In fact, as has been shown in Ref. 0, only the GDE must 
be satisfied in any case, not the Schrodinger equation. In 
the light of this fact the Weinberg program can be consid- 
ered from a new point of view: Instead of the Schrodinger 
(LS) equation, one should use the GDE. Below we show 
that the approach based on this more general equation 
allows a modification of the chiral potential model that 
can increase its predictive power. 

The leading order NN potential that has been derived 
from the Weinberg analysis of diagrams in ChPT is P| 



^^(P2,Pl) 
^.(P2,Pl)--f|| 



Kr(P2,Pl) +Cs+ CtCTi, 

\ q ■ gig • 



-Tl • T2, 



(65) 
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with q = P2 — Pi- The couphng gA is the axial cou- 
pling constant, is the pion mass, is the pion decay 
constant, and (j(t) are the Pauli matrices acting in spin 
(isospin) space. The derivation of the chiral NN poten- 
tial at higher orders was pioneered by Ordonez, Ray and 
van Kolck 0, ^| who derived a NN potential in coordi- 
nate space based upon ChPT at NNLO. Epelbaum et al. 
[l7l | constructed the first momentum-space NN poten- 
tial at NNLO. Recently the first NN potential at fourth 
order (N^LO) of ChPT was developed by Entem and 
Machleidt This potential was shown to reproduce 

the NN data below 290 MeV with the same accuracy as 
phenomenological high-precision potentials. This poten- 
tial can be written in the form 



Vsh, 



and 



(2) 



.(4) 



sh 



sh 



sh 



(66) 
(67) 



(68) 



potential in the ^Sq channel takes the form 



(P2|"^"|Pl) = Co + c2{pI+pI) + aipt+pt)- 



(69) 



At NLO one can restrict oneself to the first two terms in 
the expansion (|69|) . and the potential takes the form 



(P2|V^|pi) 



Co 



C2(P2 



(70) 



The potential ^ is very singular at short distances, 
and the LS equation with this potential makes no sense 
without regularization. Let us perform regularization by 
using a momentum cutoff. In this case the singular poten- 
tial l(7n|) is replaced by the regularized one (p2|Va|pi) = 
/*(p2/A)(p2|y|pi)/(pi/A), where the form factor f{p) 
satisfies /(O) = 1 and falls off rapidly for p > 1. At next- 
to-leading order the regularized NN potential in the ^Sq 
channel is of the form 



Vtt of the chiral potential at order t/, and yj^-* is the ly-th 
order contribution to the short-distance part. At very 
low energies the pion-exchange part of the NN interac- 
tion may be included into the contact term, and the NN 



Here is the contributions to the long-distance part (P2|V;(^^|Pi) = /*(p2/A)(^co(A)-f-C2(A)(p2+Pi))/(Pi/A)- 

(71) 

This potential has a two-term separable form and so the 
corresponding LS equation can be solved using standard 
techniques. The T-matrix obtained in this way is 



(P2|Ta(z)|pi) 



co(A) + ci(A)/4(^) + (c2(A) - ci(A)/2(z)) (pI+pI) + cl(K)h{z)plpl 



1 - (co(A) + c2(A)/4(z))/o(^) - 2c2(A)/2(z) + c2(A)/|(z) 
where the integrals In{z) are given by 

This expression can be rewritten in the form 

(P2|rA(z)|pi) = 



/*(P2/A)/(pi/A). 





C2(A)-C2(A)/2(Z)' 


{pI+pI) + cI{A)Io{z)pIpI) 


co(A) + ci(A)/4(^)' 


-1 




(l-2c2(A)/2(z)- 


fc2(A)/|(z)) 


co(A)+c2(A)/4(z)' 


-^o(^) 



r(P2/A)/(pi/A). 



I 

Taking into account that /„(z) = zm/„_i(z) -|-/„(0) and where 

'' d^~k~e^--^^\fCk)\\ (72) ^° = [co(A) +c2(A)/4(0)] [1 - (co(A) + c2(A)/4(0)) 



/2«(0) = — 



(27r)3 



and keeping only terms that contribute at leading and 
next-to-leading orders, we can write 



x/o(0) - 2c2(A)/2(0) + c2(A)/|(0)]-\ 

C2(A)fl-C2(A)/2(0)^ 



(P2|Ta(z)|pi) 

(l+C2pf)(l+C2pi)(l + 0{(g/A)n) 

Cq"^ - i^V^^Cl + 2c2zm) + J2zm' 



C2 



(73) 



co(A)+c2(A)/4(0) ' 

:72 = -2c2/o(0)-c2/2(0). 



(74) 
(75) 
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FIG. 1: The real part of the half off-shell T- matrix as a func- 
tion of the on-shell momenta p. The curves represent the real 
part of the amplitude ^(p',p) = ~{i>'\T{Ep+iO)\p) (|p'| = 10 
MeV) obtained by solving the LS equation with the potential 
1711 1 at LO (dotted line) and NLO (solid line). The exponen- 
tial regulator was employed with A = 150MeV. 
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FIG. 2: The ^^o phase shift of NN scattering. Solid dots 
represent the experimental data [l^. The dotted, solid and 
heavy solid curves represent the results obtained by using Eq. 
fM with the LEG'S C2 = 8 • 10"'^MeV"^ J2 = -0.616 and 
C4 = -2 • lO-^MeV*, J4 = -8 ■ IQ-^UeV-'^ at LO, NLO 
and NNLO, respectively. 



Note that, since the parameters co(A) and C2(A) that 
appear in Eq. (|71|l are assumed to be real, the constant 
C2 is also real. 

This is just the expression for the T-matrix shown in 
Eq. (|5U|) -with given LEC's. Thus, even if we start from 
the LS equation with the singular potential H70I) . after 
regularization and renormalization we arrive at the ex- 
pression (|50|) for the NLO 2N T-matrix we have derived 
from the requirement that the T-matrix is of the form 
(|19|l and satisfies the GDE. The constants that occur 
in the potential H7UI) should be fitted to the scattering 
length and effective range. As it follows from Eq. (|52|l . 
this means that the constants Cq is fixed by the scatter- 
ing length C'o = Ana/m, and C2 and J72 must satisfy the 
equation 

= (2c2 - ^) . (76) 

In this way, by using the regulator function of the form 

/(p/A) = exp(V/A4), (77) 

with A = 150 MeV, for C2 and J2 we get the fol- 
lowing values of the LEC's: C2 = -4.28 • IQ-^MeV-^, 
J2 = —0.463. The half off-shell T-matrix reproduced by 
the solution set fl2 of the GDE corresponding to these 
constants is shown in Figs. 1. As we see, the NLO ampli- 
tude A{p' ,p) = — (p'|T(i?p + iO)|p) is much larger than 
the leading order one even for very low momentum p'. 
This may be a signal that the EFT expansion breaks 
down in describing the off-shell behavior of the 2A'^ T- 
matrix at anomal low energies. Our finding is consis- 
tent with results obtained by Kaplan, Savage, and Wise 
__ 00 
|15|: The expansion of the NN potential = I] 

with being contact potentials breaks down at the scale 
A ~ since at this scale Vq — V2. This result is 
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FIG. 3: The real part of the half off-shell T-matrix as 
a function of the on-shell momenta p. The dotted, solid 
and heavy solid curves represent the solutions of the GDE 
for LEG'S C2 = 8 • IQ-^MeV-^, J2 = -0.616 and C4 = 
-2 ■ IQ-^UeV-'^, Ji = -8 ■ IQ-^MeV-^ at LO, NLO and 
NNLO, respectively. 



very discouraging from the EFT point of view: The am- 
plitudes should be expanded in powers of {Q/AY where 
A is a mass scale typical of the particles not included 
explicitly in the theory. In the low energy EFT, with 
the pion integrated out, one would expect A ~ m-^. The 
scale A ^ I/^/otq at which the effective expansion of 
the 2A^ potential breaks down can hardly be called a 
typical QCD scale. For this reason in Ref. it was 
proposed to expand not the potential but the Feynman 
amplitudes. This expansion was shown to converge up to 
the scale A ~ m^. However, in this way one can calculate 
only the scattering amplitudes not the off-shell T-matrix. 
The advantage of the way of the expansion of the effec- 
tive theory based on the formalism of the GQD is that 
it allows one to calculate not only the scattering ampli- 
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tudes but also the ofF-shell 2A'' T-matrix, and its EFT 
expansion converges up to the expected scale set by the 
pion mass. But for this we must not restrict ourselves to 
the solution sets associated with potential model when 
the parameters C2 and J2 are constrained not only by 
Eq. but also by Eq. (jJH). Indeed, in general these 
LEC's determining the solution set ^2 must satisfy only 
Eq. H77|) . This means that fitting to the scattering data 
does not allow one to obtain these constants. However, 
it is natural to expect that the off-shell T-matrix must 



approach order by order to the true one with the same 
accuracy as the on-shell matrix elements that determine 
the phase shift. This puts some additional constraint on 
the values of the constants. Let us investigate the prob- 
lem of the convergence of the effective expansion up to 
NNLO. At this order the 2iV T-matrix is given by Eq. 
(|56|l . For simplicity we will restrict ourselves to the case 
when the function <i>o(P2, Pi) in Eq. (|64|) is equal to zero, 
and the expression for the NNLO 2N T-matrix in the ^Sq 
channel takes the form 



(P2|tW(z)|Pi) 



(1 



C2VI 



C2P1 



■ C4p1) 



C^^ - ^V^i^Hfl + 2Rec2zm 



2Rec4 



x(l + 0{QVAn), |Pi|<A, |P2|<A. 



(78) 



This is a good approximation to the 2A^T-matrix because 
it has a pole near £^ = in the ^5'o channel. The re- 
sults of our calculations have shown that an acceptable 
convergence of the effective expansion is achieved in the 
case when the parameter C2 is positive. For example, 
as we see from Figs. 3 and 4, convergence up to the ex- 
pected scale takes place for the following set of the LEC's: 
Co = -0.0016MeV-2 (fixed at LO); ca = 8 • IQ-^MeV-^ 
and J2 = -0.616 (fixed at NLO); C4 = -2 • lQ-^MeV~^ 
and J4 = -8 • IQ-^Me^^"^ (fixed at NNLO). Note that 
in this case we deal with the set of the solutions of the 
GDE that cannot be associated with any potential. This 
is because the condition H75(l that can be satisfied in the 
case of the solutions obtained starting from the LS equa- 
tion with the singular potential H7U|) allows only negative 
C2. In fact, taking into account that in the ^6*0 chan- 
nel iaro = -(35M^)^ while C2 is of order ^ ~ i^f^ 
from Eq. H76(l we conclude that J2 is close to — iaroCj^^, 
and hence is negative. From Eqs. H72|l and H75|) than it 
follows that the constant must be negative. 

For momentum p ^ TOtt explicit pion fields must be 
included in the theory. In this case in analyzing the time 
ordered diagrams for 27V T-matrix one has to take into 
account that the contributions from the " 2 irreducible" 
diagrams manifest themselves not only as contact interac- 
tions but also as pion-exchange ones. In this context the 
2A^ T-matrix obtained in the previous section should be 
considered as a contribution from the processes in which 
only contact interactions come into play. According to 
the Weinberg proposal this contribution should be ob- 
tained by iterating the corresponding contact potential 
in the LS equation. On the other hand, the solutions of 
the LS equation have the following asymptotic behavior 

(p2|T(z)|pi) -> (P2|F|P1), 

\z\ — ^oo 

which means that in the limit of high energies the T- 



matrix gets the dominant contribution from the poten- 
tial. In our case, when this T-matrix describes the sum 
of contributions from the all relevant contact diagrams in 
ChPT and the potential is its "2A^ irreducible" part, the 
above should mean that this part gives the main contri- 
bution in the limit of high energies. However, this is not 
the case because all Feynman diagrams with Cq's at each 
vertex give a contribution of the same order. In the chiral 
potential approach this problem is solved by putting in 
ad hoc short-ranged form factors in front of the contact 
interactions. As we have demonstrated, in the ^S'o chan- 
nel at NLO after renormalization of the solution of the 
LS equation with potentials regularized in this way we 
arrive at the expression (|49|l we have obtained, by using 
the GDE and the momentum dependence of the above 
sum of the contact diagrams that is given by Eq. I|19() 
without extracting any potential. It can also be shown 
that renormalization of the solution of the LS equation 
with the NNLO contact potential gives rise to the expres- 
sion (|64l) . However, as we have seen, there are solutions 
of the GDE consistent with symmetries of QCD that can- 
not be obtained by starting from the LS equation with 
regularized contact potentials. This is not surprising for 
the following reasons. 

As is well known, the short-distance physics in ChPT 
is not the full short-distance physics of QCD. At dis- 
tance scales ~ 1/A^ mechanisms not encoded in ChPT 
will come into play, and will regulate the behavior of the 
nuclear forces at short distances. On the other hand, in 
quantum field theory, if we spread interaction in space, 
we spread it in time as well. Thus the processes that in 
describing low energy nucleon dynamics should be used 
as a short-range part of the NN interaction are nonlocal 
in time. In the ordinary quantum mechanics nonlocality 
in time of an interaction, when it is used as a "funda- 
mental" interaction generating the dynamics of a system 
creates the conceptual problems with causality and uni- 
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tarity. This is because the Shrodinger equation is local 
in time, and the interaction Hamiltonian describes an in- 
stantaneous interaction. As we have seen, the formalism 
of the GQD provides the extension of the quantum dy- 
namics needed for solving this problem, and the GDE 
allows one to describe the evolution of a system with a 
nonlocal-in-time interaction. It is remarkable that such 
a nonlocalization can take place if and only if the UV be- 
havior of the matrix elements of the evolution operator 
is "bad". 

As we have shown, the requirement that the 2N T- 
matrix satisfies the GDE and has momentum behavior 
shown in Eq. (|19|l allows one to construct it order by 
order in the effective expansion. Note that at any or- 
der (p2|T(^-'^^(z)|pi) constructed in this way satisfies the 
GDE and hence is unitary. The corresponding interac- 
tion operator describes the effective interaction that is 
spread not only in space but also in time. This man- 
ifests itself in the momentum and energy dependence 
of {'P2\B^eff\z)\pi) ■ The regularization of the contact 
potentials spreads the interaction described by this po- 
tential only in space. Iterating this potential in the LS 
equation spreads the interaction also in time. So in the 
energy region relevant for describing low energy nucleon 
dynamics the effective interaction reproduced by the reg- 
ularized potential is nonlocal in time. However, since this 
potential gives the main part of the contribution to the T- 
matrix only in the region much above the scale A^ where 
the nonrelativistic two-nuclcon LS equation is physically 
meaningless, in this case the solution of the LS equation is 
only of formal importance. To clarify this point note that 
the LS equation is equivalent to the Schrodinger equation 
that implies the existence of the infinitesimal generator 
of the group of the evolution operator V{t). However, for 
the operator {y{t,Q) — l)/t to approach to the limiting 
one sufficiently close, one has to cross the region of times 
much below the time scale 1 /A^. For such times, as it fol- 
lows from Eqs. (|13|l and (|14|l . the dominant contribution 
comes from the operators T(z) corresponding to z much 
above the scale A^. This means that the dynamics is 
not governed by the Schrodinger equation formulated in 
the terms of low energy degrees of freedom. The above 
might be considered as an evidence of inconsistency of 
the cutoff approach, but the formalism of GQD provides 
a new reason to use the LS equation with a regularized 
potentiahThe solution of the LS equation with any reg- 
ularized potential satisfies the condition H59() . and hence 
the solution for z e I? can be used as the effective in- 
teraction operator in the GDE. In this case the existence 
of the infinitesimal generator of the group of the evolu- 
tion operators is not required. However, as we have seen, 
not all the solutions of the GDE relevant for the prob- 
lem can be obtained in the above way. In other words 
there are the interaction operators satisfying the condi- 
tion (|59|l that cannot be represented as a formal solution 
of the LS equation with any regularized potential. The 
above means that it is possible that not all the effects 
of short-distance physics of QCD can be incorporated 



in the effective theory of nuclear forces in the standard 
way where the short-range part of the NN interaction 
is described by the regularized contact potentials. This 
motivates a modification of the chiral potential model 
that could incorporate the description of the short-range 
interaction based on the use of the GDE. 

In the theory with pious the interaction operator 
should contain terms describing both the short-ran ge a nd 
long-range parts of the NN interaction. In Ref. [l^ it 
has been shown that at leading order in the ^5*0 chan- 
nel the NN interaction operator that is determined by 
the behavior of the 2N T-matrix in the energy region V 
where the dominant contributions to the 2N T-matrix 
come from the "2N irreducible" diagrams of ChPT is of 
the form 

(P2|i?(°)(z)|pi) = (P2|i?l°^(z)|pi) + (P2|T^i'')|Pl), (79) 

where the potential V^i^-* is given by Eq. (jSHJ), i.e., is just 
the LO component of the long-range part of the chiral 

potentials, and (z) describes the short-range compo- 
nent of the NN interaction at LO and is given by Eq. 
(|33|l . In the present paper we cannot discuss the theory 
with pions in detail. However, the arguments presented 
in Ref. lead us to the conclusion that at any order 
the 2N interaction operator can be represented as a sum 
of the nonlocal-in-time short-range component being the 
sum of the relevant contact diagrams and instantaneous 
part that is the sum of the irreducible pion-exchange di- 
agrams 

(P2|i?(z)|pi) = (P2|i?.h(z)|pi) + (p2|y.|pi), (80) 

where the long-range component (p2|K-|Pi) uiay be de- 
scribe by that of the chiral potential at the corresponding 
order. The solution of Eq. (jSJ with this interaction oper- 
ator is not more complicated than the solution of the LS 
equation with an ordinary potential. In fact, the solution 
of Eq. (|SJ| with the boundary condition @ where B(z) 
is given by Eq. (|80|l can be represented in the form 

T{z) = Tsh{z) + (l + T,,,(z)Go(z))f,(z) 

X {l + Ga{z)Tshiz)), (81) 

where Tsh{z) is the T-matrix describing the dynamics by 
Bsh{z) alone, i.e., is the solution of equation 

dTshjz) _ f (fk T,fe(z)|k)(k|T,^(z) 

dz J (27r)3 {z~Ek+iO)^ ' ^ ' 

with the boundary condition Tsh{z) Bsh{z), and 

\z\ — * oo 

Ttj{z) satisfies the equation 

f^{z) = V^ + V^Gsh(z)f^{z), 
Gsh (^) = Go (z) + Go {z)T,h (^)Go (z) . (83) 

Note that, by using the "two-potential trick" the so- 
lution of the LS equation with the potential V = Vsh + V-n 
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can be written just in the same form. The only difference 
is that in the latter case Tgh{z) is the solution of the LS 
equation with the potential Vsh{z), while in the general 
case where Eq. (|83ll is derived from the GDE with the 
interaction operator H80() . Tsh{z) is the solution of Eq. 
(|83|l with the interaction operator Bsh{z)- 

At NNLO the NN interaction in the ^ 5*0 channel can 
be parametrized by the operator 

(P2|i?(z)|pi) - (P2|i3g^^(z)|pi) + (P2|^^.|P1>, (84) 

(P2|^.|P1) = (P2|K^°)|Pl) 

+ (P2|yi2)|pi) + (p2|yi^)|pi), (85) 

where the operator b'^j}j:{z) is given by Eq. H6U|) . and 

(P2|^7r''''|Pi) are the pion-exchange components of the 
chiral potential (|66|l projected to the ^5*0 state. Thus for 
the NNLO 2N T-niatrix in the ^Sq channel we have 

(p2|T(z)|pi) = (P2I {T^l\z) + T^^{z)Goiz))f^iz) 

X (1 + Go{z)T^t\z)) Ipi) (1 + {{Q/Af]) , (86) 

where T^iJ{z) is given by Eq. H49|) . At the same time, 
as we have shown, the result of resuming the contact 
potential (|70|l to all orders can be represented in the form 
(innj) with the LEG'S given by Eqs. GH) and lO. In other 
words, there are the solution sets at which one can arrive 
starting from the LS equation with the contact potential 
(|70|l . and in the case when T^1\z) in Eq. (|85|l belongs 
to such sets this equation represents the T-matrix (more 
precisely a set of the T-matrices) that coincide with the 
2N T-matrix obtained by solving the LS equation with 
the NNLO chiral potential. 

Thus in the particular case the GDE with the inter- 
action operator can reproduce the same results as 
the LS equation with the NNLO chiral potential. How- 
ever, Eq. (|86|) that represents the solutions of the GDE 
with the interaction operator (|84|l allows a much wider 
class of the solutions, and only experiment can discrimi- 
nate between them. As we show in the next section, the 
above freedom in choosing values of the LEC's character- 
izing the short-range part of the A^'A'^ interaction allows 
one to significantly change the off-shell predictions of the 
model and at the same time to keep the A^A^ phase shifts 
unchanged. This gives us the hope that the above modi- 
fication of the chiral potential model could increase their 
predictive power in describing few-body and microscopic 
nuclear structure problems where the off-shell properties 
of the A^ A^ interaction play an important role. 

VI. EFFECTS OF THE SHORT-RANGE 
FORCES ON THE OFF-SHELL 2N T-MATRIX 
AND THE PP BREMSSTRAHLUNG 

One of the simplest processes involving the half off- 
shell A'^A'^ interaction is proton-proton bremsstrahlung. 



Many years ago, it was suggested to use the 
bremsstrahlung process as a tool to discriminate between 
the various existing 2N potential models However, 
it has been shown that once all ingredients which are 
necessary for the calculation of the bremsstrahlung pro- 
cess are introduced, the predictions of various modern 
potential models do not differ significantly, independent 
of whether they are finite range or energy dependent 
[2^ . The problem is that there is a significant discrep- 
ancy between these predictions and experimental data of 
KVI I3 for asymptotic proton angles: As was shown in 
Ref.|8j, the microscopic calculations reproduce the pp'^ 
cross sections reasonable well, however, for certain kine- 
matics for which the final pp system has a low kinetic 
energy, there is a significant over predictions of the data 
by about 20% — 30%. Since the discrepancy between the 
results of the microscopic calculations and the experi- 
mental results increases as the kinetic energy of outgo- 
ing protons decreases, one would speculate that including 
the Coulomb force might resolve the problem. However, 
it has been shown by Cozma et al. [3 that this is not 
the case. An important part of this discrepancy between 
theory and experiment originates in a poor description 
of the NN interaction in the -^Sq channel at low energies 
0. However, the above refers only to potentials model, 
while, as we have seen, there are such short-distance ef- 
fects that can be represented by the GDE but not by 
any potential model. This leads us to the conclusion 
that the possible source of the discrepancy between the 
data and the microscopic calculations is the inability of 
the potential approach to incorporate such effects. Let 
us now examine this possibility in the same way that was 
used by Gozma et all for investigating the effect of in- 
cluding the Goulomb interaction. In Ref. 9] this effect 
was studied by using the separable approximation to the 
A^A^ interaction. The form factor ^(p) in the separable 
potential 

^(P2,Pi)-A5*(p2).g(pi) (87) 
was chosen in the form 

5(P) = 2 I /?2 ' 
-I- 

and the parameters A and (3 were fixed by the scatter- 
ing length a and the effective range tq. The Goulomb 
effect was investigated by adding the Goulomb poten- 
tial to the separable potential, and these corrections was 
shown to be minor in the region needed for the KVI 
bremsstrahlung experiment. 

Let us now, instead of the Goulomb effect, investigate 
the effects of the including the short-range NLO interac- 
tion operator i?^^ (z) given by Eq. (|60|l in combination 
with the potential (|87|l . In this case the interaction op- 
erator of our problem is of the form 

(p2|i?e//(z)|Pl) = (P2|Sf,^(z)|pi) 

+Ag*(P2)g(Pi). (88) 
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For such a choice of the interaction operator, Eq. H83() leads to the following T-matrix: 



(P2|T(z)|pi> 
with 



tsh{z) = 

t^{z) = 

IiTshiz) = 



1 



l~~tsh{z)U{z)I^sh{z)i;,,,{z*) 



(l + 0{QVAn) 



(89) 



Cq' - ^./^^[l + 2Rec2zm] + J2zm 
1 

T^^i r d-'fc |3(k)~ 



(Pk g(k)V'o(k) r m ^~zm - /3 



(27r)3 z-£:fe L47r 1 + ZTO//32 

I 



1 + cXzm] — Jj 



l + OiQ^A'}), 



(90) 



where 

Xsh = ^ y',g(k)(v*(k) - (91) 

Here we restrict ourselves to the case where the function 
■00 (p) that occur in Eq. H49|l is real. Equation l|89|) rep- 
resents a set of the solutions of the GDE that coincide 
with the true 2N T-niatrix with an accuracy of order 

0{Q^ /A^}, because, as we noted, the effective operator 
(2) 

^sh i-^) represents a set of the operators describing the 
short-range part of the NN interaction. Note that this 
solution set is determined not only by the LEC's Co, C2, 
J^2, /?, and A that occur in the interaction operator H88() 
but also by an additional constant J^-^sh- This means 
that really the operator H88|l being the sum of the short- 
range interaction operator and the potential V^r is not so 
close to the true solution in the domain T) to uniquely 
determine the set of the solutions of the GDE that co- 
incides with an accuracy of order 0{Q^ / h^}. In Ref. 
[19] this situation was explained by using the example of 
the interaction operator H79|l generating low energy nu- 
cleon dynamics at leading order. As it has been shown 
in Ref. 19], at high z this interaction operator is not 
so close to the true 2N T-matrix to uniquely determine 
the relevant solution. This means that the operator H79() 
does not satisfy the condition (|59|) . By substituting the 



operator (|79|l into Eq. (|59|l . one can obtain terms that 
should be added to the short-range part of the interaction 
operator In principle the interaction operator H88() 
should be corrected in the same way. However, one may 
treat the interaction operator of the form H88|l keeping 
in mind that this operator does not determine a unique 
solution of the GDE (more precisely a unique solution 
set). Indeed, Eq. (|89(l represents a family of the solu- 
tion sets: In general the solutions belonging to the sets 
with different values of the parameter J7irs/i given by Eq. 
(|91ll do not coincide with each other with an accuracy 
of order 0{Q^ / K^} . So to select one of these sets we 
should chose the parameter Jt^sK- In general this param- 
eter should appear in the above hidden additional terms, 
and the corrected interaction operator should determine 
a unique solution (solution set) of the GDE. Thus the 
set ^2 of the solutions of the GDE that coincide with the 
true 2iV T-matrix with an accuracy of order 0{(3^/A'^} is 
determined by the LEC's Cq, C2, J2, J-ksH characterizing 
the short-range part of the NN interaction and the con- 
stant /3, and A parametrizing its long-range component. 
Each of the solutions belonging to this set corresponds 
to some definite function V'o(p) satisfying the conditions 
(|77|l . llSni) and linU with given LEC's C2, J2 and J-r^sh, 
and for momenta below A can be represented as 



(p2|T(z)|pi) = 



(1 -I- cIpI){1 + C2pl)tsh{z) + g*(P2)g(Pl)t^(^) 
1 - tsh{z)U{z)Il^f^{z*)I^,h{z)) 

(1 + c5pi)g(pi)/:,Jz*) +5*(P2)(1 + C2P?)/..,.(z))t.^(z)t,(z) 



tsh{z)t^{z)I*Az*)hsh{z) 

I 



i + o{gVA=^} 



(92) 



The constants of the model can be obtained from the re- 



quirement that Eq. (|92|l reproduces observables with an 
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accuracy of order 0{(5'^/A^}. Obviously the fitting to 
the pp scattering data is not sufficient for this, because 
these constants determine not only the on-shell but also 
the off-shell properties of the 2N T-matrix. For their de- 
termination the comparison with experimental data such 
as the bremsstrahlung data in which the off-shell proper- 
ties of the NN interaction manifest themselves is needed. 
For example, one can use the T-niatrix H92|) for describing 
the elastic 2N T-matrix in the pp'y model of Martinus et 
al. [23, 01 ) and try to fit the constants to reproduce the 
experimental cross sections. However, in this paper we 
restrict ourselves to the investigation of the effects of the 
short-range interaction on the off-shell 2N T-matrix in 
the region needed for KVI bremsstrahlung experiment. 

In Figs. 4,5,6, and 7 the results of the calculations 
of the T-matrix given by Eq. (|92|l . The parameters A 
and (3 characterizing the long range component of the 
interaction operator H88|l were determined by fitting the 
scattering amplitude reproduced by the pure potential 
model to the scattering length a = — IS.lfm and the ef- 
fective range ro = 2.60fm. These values of a and tq 
correspond to the case when the Coulomb interaction is 
switched off Q. The results of our calculations plotted 
in Figs. 4 and 5 show that for some values of the LEC's 
C2, J2, and JTrsh the difference between the half off-shell 
amplitudes {pi\T {Ep + iO)\p) with and without the short 
range corrections becomes larger as the momentum go- 
ing from the on-shcU point, while, as we see from Fig. 
4, including the contact term keeps the phase shift un- 
changed. In contrast, as has been shown in Ref. [sl, 
the Coulomb corrections become minor for high off-shell 
momenta, and hence cannot remove the discrepancy be- 
tween theory and experiment in describing the pp^ pro- 
cess. Our results show that with an appropriate choice 
of the constants the short range part of the NN interac- 
tion gives a significant contribution to the 2N T-matrix 
in the case when the laboratory kinetic energy i?iab of 
the outgoing protons is below 15 MeV and the off-shell 
momentum is sufficiently high. This is just the region 
needed for the KVI bremsstrahlung experiment. A 
large discrepancy between experiment and the predic- 
tions of the potential models appears in the kinematical 
regions where the cross section has peaks. The presence 
of these peaks is the result of strong final-state interaction 
at small values of the kinetic energy of outgoing protons 
(-Blab < 15MeV). This means that the over prediction of 
the potential models in describing the ppj cross sections 
may be a reflection of the fact that really the NN inter- 
action at low energies is not so strong as it is predicted 
by the potential models. In Figs. 6 and 7 we plot the 
relative difference 



A(^lab,p') = 



{p'\TAEj, + iO)\p)\-\{p'\TiEp + iO)\p) 
\{p'\T{Ep + tO)\p)\ 



Eiah — 2p^ /m, between predictions of the model with 
and without the short-range component for the half off- 
shell T-matrix where T^{z) is the solution of the LS 
equation with the potential H87|) . As we see, for high 
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FIG. 4: A comparison between the pp phase shifts (^5*0 chan- 
nel) for the potential 18711 fitted to a = —18.1 fm and re ~ 2.60 
fm (dotted line) and the predictions of the model corrected 
by including the short effects (solid line). 
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FIG. 5: The T-matrix as a function of the off-shell momentum 
p at a laboratory energy of 10 MeV. The solid and dashed 
curves represent, respectively, the real and imaginary parts 
of the normalized T-matrix {p\Tn{Ek + iO)|fc)/|(0!r^(£;fc -I- 
iO)|fc)|, where T-jriz) is the solution of the LS equation with 
the potential I87I I. The heavy solid and bold type dashed 
curves represent, respectively, the real and imaginary part of 
{p\T{Ek + iQ)\k) /\{Q\T^{Ek + iO)|fc>|, where T{z) is given by 
Eq. The parameters are following (5 = 241.468MeV, 



A 



-1.0269 ■ 10-*MeV-^ Co = 1.29 • 10-''MeV-^ X 



-17942MeV^ 



C2 



2.3 ■ 10-^MeV-^ J2 = 2.75. 



off-shell momenta, the magnitude of A(i!^iab,p) charac- 
terizing the effect of the short-range forces on the 2N 
T-matrix grows rapidly as the energy _Eiab going to zero 
in the region Euh < 15 MeV. Note that in this kine- 
matic region the difference between the predictions of 
the microscopic models for bremsstrahlung and the KVI 
experimental data has just the same energy dependence 
8J. This leads us to the conclusion that the improve- 
ment of the low-energy part of the strong interaction that 
can be achieved by correcting its short-range component 
may remove the discrepancy between theory and exper- 
iment. However, for including the short-range compo- 
nent to give rise to the above effect, the constants C2 and 
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FIG. 6: The difference A(i?iab,p) between the predictions of 
the model with and without the short-range corrections for 
the T-matrix in the region needed for the KVI experiment 
(p = 300 MeV). 
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FIG. 7; Same as Fig. 6 but for p = 500 MeV. 



J^2 should have such vahies that are inconsistent with 
any potential model. In fact, if the LEC's C2 and J2 
correspond to some contact potential then they must 
be related by Eq. (|75|l . From this it follows that for 
C2 = 2.3- 10-'^MeV"2 ^nd J2 = 2.75 which lead to the re- 
sults shown in Figs. (4-7) the integral / dk\f{k)\'^ should 
be of order - = 1.26 • 10^A-\ where the cutoff is 

set by the chiral symmetry breaking scale A^^, — IGeV. 
However, such a size of this integral is inconsistent with 
the fact that the regulator function /(fc/A) must satisfy 
/(O) = 1 and fall off rapidly for k/A > 1. This means 
that the effective 2N T-matrix (|92|l corresponding to the 
above LEC's cannot be associated with any potential. 

The above model can be considered as a combination 
of the standard nuclear physics approach based on the 
use of the NN potentials and the effective theory of nu- 
clear forces. As we have seen, combining our approach to 
the description of the short-range part of the NN inter- 
action with the separable potential model yields very en- 
couraging results in investigating possible origins of the 
discrepancy between the theoretical predictions for the 



pp bremsstrahlung and experiment. In the same way one 
can construct a more realistic model by replacing the po- 
tential (|87|l in Eq. H88|) with one of the high precision 
phenomenological NN potentials [25-30]. These poten- 
tials contain the terms describing the short-range part of 
the NN interaction, and hence in this case the operator 
Bsh{z) in Eq. (|79|l should be considered as a correction 
to the model dependent short-range components of the 
realistic potentials in which the ambiguous short-distance 
parametrization is used. It is important that such cor- 
rections may allows one not only to remove or minimize 
the model dependence but also improve the description 
of the off-shell properties of the NN interaction that is 
needed for removing the discrepancy between the theory 
and the experiment in describing the pp bremsstrahlung. 



VII. THE OFF-SHELL BEHAVIOR OF THE 2N 
r-MATRIX AND THE THREE-NUCLEON 
PROBLEM 

As we have shown, our formulation of the effective the- 
ory of nuclear forces allows one to construct, as an in- 
evitable consequence of the basic principles of quantum 
mechanics and symmetries of QCD, not only the 2N scat- 
tering amplitude (in this case we reproduce all results of 
the standard subtractive EFT approach), but also the 
off-shell T-matrix, and hence the evolution and Green 
operators. The advantage of this formulation over the 
tradition potential approach is that it allows one to take 
into account the constraints on the off-shell properties of 
the NN interaction put by the symmetries of QCD. This 
is very important, because, as is well known, the off-shell 
behavior of the 2A^ T-matrix may play a crucial role in 
solving the many-nucleon problem and is an important 
factor in calculating in-medium observables j3lj | and in 
microscopic nuclear structure calculations. This results, 
for example, in the fact that the predictions by the Bonn 
potential for nuclear structure problems differ in a char- 
acteristic way from the ones obtained with local realistic 
potentials . The off-shell ambiguities of realistic NN 
potentials are argued to be one of the main causes of 
many problems in describing three-nucleon systems. 

The modern potential models that describe 2N scat- 
tering data to high precision can not guarantee that a 
similar precision will be achieved in the description of 
larger nuclear systems. In fact, the simplest observable 
in the 37V system, the binding energy of the triton, is un- 
der predicted by the realistic NN potentials which are so 
successful in describing the 27V observables. The energy 
deficit ranges from 0.5 to 0.9 MeV and depends on the 
off-shell and short-range parametrization of the 2N force 
[33! . In order to resolve this problem one has to take into 
account three-nucleon force (3NF) contributions to the 
3iV binding energy. The common way of solving the 3A'^ 
bound state problem is to use in the Schrodinger equation 
phenomenological NN potentials and then to introduce 
a 3NF to provide supplementary binding. However, from 
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the point of view of the three-nucleon problem, it is not 
sufficient to generate a phenomenological NN potential 
that perfectly reproduces the 2N scattering amplitudes. 
One must also generate a NN potential by using theo- 
retical insight as much as possible in order to constrain 
the off-shell properties of the 2N T-matrix. If this is 
not the NN potential which fits precisely the 2N 

phase shifts but produces the erroneous off-shell behavior 
of the T-matrix would not provide reliable results for the 
3A^ system, nor can be used to test for the presence of 
3 A'' forces. 

The models for the 3NF that are usually used for solv- 
ing the problem with the triton understanding are based 
on two-pion exchange with intermediate A-isobar excita- 
tion. However, these 3NF models cannot explain the Ay 
puzzle which refers to the inability to explain the nucleon 
vector analyzing power Ay in elastic nucleon-deuteron 
(Nd ) scattering at low energy. As has been shown in Rcf. 
|33 | , it is not possible with reasonable changes in realistic 
NN potentials to increase the Nd Ay and at the same 
time to keep the 2N observables unchanged. The same 
situation also takes place in the case of chiral potentials 
|35j |. This means that the introduction of three-nucleon 
forces (3NF) is needed for resolving the problem. How- 
ever, as it turned out, conventional 37VF's change the 
predictions for Nd Ay only slightly and do not improve 
them [36-38]. This motivated new types of 37VF's [39- 
41]. On the other hand, reliable 3A^ calculations and 
even testing for the presence of 37V forces require to con- 
strain the off-shell properties of the 2N. Moreover, these 
properties play a crucial role in a new 3A'^F proposed by 
Canton and Schadow In Ref. [I^ this 3A^F has 

been suggested as a possible candidate to explain the Ay 
puzzle. The 3iVF contribution from this force is fixed 
by the 27V T-matrix describing the underlying 2N inter- 
action while the pion is " in flight" . The expression for 
this 3A^F derived in Ref. ^ contains the off-shell 2iV 
T-matrix, more precisely its subtracted part 

?12(P2,Pl,2;) = il2(P2,Pl,z) - 'yi2(P2,Pl), (93) 

where pi and P2 are Jacobi momenta of nucleons 1 and 
2, and the potential - like term Vi2(p2, Pi) contains only 
OPE/OBE - type diagrams. The subtraction in Eq. 
is needed to take into account a cancellation effect which 
has been observed in Refs. 0, This cancellation 
involves meson retardation effects of the iterated Born 
term, and the irreducible diagrams generated by sub- 
summing all time ordering diagrams describing the com- 
bined exchange of two mesons amongst the three nucle- 
ons. In principle there are no free parameters to adjust, 
and the 3A^ force is completely determined by the 2N 
T-matrix. However, as has been shown in Ref. [4^ . in 
order to explain the Ay puzzle, instead of the subtracted 
T-matrix shown in Eq. H93|l . one has to use the amplitude 
defined according to the prescription 

il2(P2,Pl,2) = c(z)ti2(p2,Pl,2) " 'yi2(P2,Pl), 

with the effective parameter c(z), which represents an 



overall correction factor for the far-off-shell 2N T-matrix. 
Ideally, this parameter should be one for the 2N potential 
to provide a reliable extrapolation of the 2N T-matrix 
down to z ~ —160 MeV. However, as has been shown 
in Ref. ^3l> none of the existing 2N T-matrices can 
guarantee the off-shell behavior that is needed for the 
explanation of the Ay puzzle with the parameter c{z) set 
to one, since they are all constrained by such data the 
deuteron pole and the 2N phase shifts. For example, in 
order to reproduce the nd experimental data with the 
Bonn B potential, the factor c{z) must be set to 0.73 for 
the energy 10 MeV. 

The above gives us the expect that the origin of dis- 
crepancy between the theory and experiment in describ- 
ing Ay is the same as in the case of the pp bremsstrahlung: 
The potential models cannot reproduce the relevant off- 
shell behavior of the 2N T-matrix. This should mean 
that the answer lies beyond these models. The advan- 
tage of the formulation of the effective theory is that it 
allows one to constrain the possible T-matrix by the sym- 
metries of QCD, but at the same time, permits the solu- 
tions that cannot be reproduced by any potential model. 
This opens new possibilities to resolve the Ay puzzle. 
Note that the generalization of the above results to the 
P waves needed for solving this problem is straightfor- 
ward. 



VIII. SUMMARY AND DISCUSSION 

We have shown that from the Weinberg analysis of 
time-ordered diagrams for the 27V T-matrix in ChPT it 
follows that nucleon dynamics at low energies is governed 
by the GDE with a nonlocal-in-time interaction operator. 
A remarkable feature of the GDE which follows straight- 
forwardly from the first principles of quantum mechanics 
is that it allows one to construct all physical amplitudes 
relevant for the theory under consideration by using the 
amplitudes describing processes in which the duration 
time of interaction is infinitesimal. It is natural to as- 
sume that the most of contribution to these amplitudes 
comes from the processes associated with a fundamental 
interaction in a quantum system. This point manifests 
itself in the boundary condition Q. If we do not con- 
sider a theory that is valid up to infinitely high energies 
(infinitesimal times), then the duration times of interac- 
tion much below the time scale of low energy physics but 
much above the scale of underlying high energy physics 
should be considered as " infinitesimal" , and the opera- 
tor describing the processes with such duration times of 
interaction can be used as an interaction operator. This 
in turn means that the amplitudes describing this "fun- 
damental" interaction in the low energy theory can be 
computed in terms of the underlying high energy physics. 

The GDE can be represented in the form of the differ- 
ential equation ||HJ) for the operator T{z). The boundary 
condition on this equation is shown in Eq. (Q, where 
the operator B{z) describes the fundamental interaction 
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in this system. By definition, this operator must be so 
close to the true T-matrix in the hmit \z\ — > oo that 
the GDE with the boundary condition have a unique 
solution. The above means that really this region of " infi- 
nite" energies with which we have to deal is the domain T) 
that lies much above the scale of the low energy dynamics 
but much below the scale of the underlying high energy 
physics. Correspondingly for describing the low energy 
dynamics we have to start with the boundary condition 
(|58|l that implies that the interaction operator B{z) is 
so close to the relevant T-matrix inside the domain T) 
that Eq. (jS)) with this initial condition has a unique so- 
lution. In the theory of nuclear forces this domain is 
a region of energies that are high enough for the most 
of contribution to the 2N T-matrix to come from pro- 
cesses that are described by the irreducible 2N diagrams 
for the 2N T-matrix (these processes are associated with 
a "fundamental" interaction), but not so high for the 
heavy degrees of freedom manifest themselves explicitly. 
By using the analysis of the time-ordered diagrams for 
the 2N T-matrix in ChPT inside the domain V where 
the structure of the theory is much simpler than in the 
low energy region, we can obtain the NN interaction op- 
erator which should be used in the boundary condition 
(I58|l on Eq. JH)) . We have shown that the dynamics which 
is generated by the interaction operator obtained in this 
way is non-Hamiltonian. This is because, for low energy 
nucleon dynamics to be Hamiltonian, the operator T[z) 
must have a negligible dependence on z inside the domain 
T), while, as we have shown, this is not the case. In other 
words, in the nonrelativistic limit QCD leads through 
ChPT to low energy nucleon dynamics that is not gov- 
erned by the Schrodinger equation. However, this does 
not mean that the low energy predictions of QCD are 
not consistent with quantum mechanics: Only the GDE 
must be satisfied in all the cases. We have shown that 
the requirement that the 2N T-matrix satisfies the GDE 
as a well defined equation free from UV divergences and 
has the structure consistent with the symmetries of QCD 
allows one to construct it order by order in the effective 
theory expansion without resorting to regularization and 
renormalization. 

The above point is illustrated in detail by using the pi- 
onless theory of nuclear forces. For the ^Sq channel 2N 
T-matrix that, as it follows from the Weinberg analy- 
sis, should be of the form (|19|l . we have derived Eq. (|23|l . 
This equation allows one to obtain the 2N T-matrix with 
an accuracy of order 0{(Q/A)^'^+^} if it is known with 
an accuracy of order 0{{Q / h.Y-^~^} . In this way at any 
order one can obtain the most general possible 2N T- 
matrix consistent with the symmetries of QCD. At NLO 
this T-matrix is given by Eq. (|49|l . This equation repre- 
sents the set VI2 of the solutions of the GDE that coincide 
with the true 2A'^ T-matrix with the accuracy up to order 
0{Q^ / A?}. This solution set is generated by the effec- 
tive interaction operator (|60|l . The fact that we have 
to deal with the set of the solutions means that at this 
order our knowledge about the details of the NN inter- 



action is not sufficient for discriminating between them. 
The set ^2 is determined by the EEC's Co, C2, and J2 
that parametrize the effects of the high energy physics 
on low energy nucleon dynamics. It is shown that, the 
2N T-matrix obtained in this way leads to the KSW ex- 
pansion of the scattering amplitude 14] shown in Eq. 
(|57|l . However, in contrast with the standard subtrac- 
tive approach to the EFT of nuclear forces where this 
expansion is obtained by summing the bubble diagrams 
and performing regularization and renormalization, we 
reproduce the same results in a consistent way free from 
UV divergences. At NNLO the 2A^ T-matrix is given 
by Eq. H54|l . This equation represents the set $^4 of 
the solutions of the GDE that coincide with the true 2N 
T-matrix with an accuracy of order 0{Q^ / h.^} . In the 
same way one can obtain the solution set U,2j\f that co- 
incide with the true 2N T-matrix with the accuracy up 
to order 0{{Q/Kf^+^}. As M increases the set VL2j\r 
becomes smaller and smaller, i.e., one approaches to the 
true 2A^ T-matrix closer and closer. 

The above strategy can be also used for constructing 
the T-matrix in the theory with pions. Another way of 
applying the results obtained in the paper is to use them 
in combination with the chiral potential approach. We 
have shown that at NNLO the NN interaction in the 
^5*0 channel can be parameterized by the operator H84|l . 
The solution of the GDE with this interaction operator 
is given by Eq. In the particular case when the 

short-range part of the 2A^ T-matrix T),^' {z) that enters 
in Eq. (|86|l can be associated with some potential model 
this equation leads to the same results as the chiral model 
with the NNLO potential. At the same time, as we have 
seen, there are such EEC's characterizing the short-range 
part of the interaction operator (|84|) in the case of which 

Tgf^ {z) cannot be associated with any potential model, 
and Eq. H86|l gives rise to the 2A^ T-matrix that cannot 
be reproduced by the LS equation with the NNLO chiral 
potential. Thus, being equivalent to the chiral potential 
model in a particular case, above model allows one to 
go beyond the potential models, and such an extension 
of the standard chiral approach may be necessary to re- 
move the discrepancy between the theory and experiment 
still existing in nuclear physics. This point has been il- 
lustrated by using the example of the model where the 
long-range part of the interaction is described by the po- 
tential H87|l . This is the model that in Ref. 9] was used 
for investigating the effects of including the Coulomb in- 
teraction on the pjTf process where the Coulomb poten- 
tial is replaced by the short-range interaction operator 
(|88|l . By using this model we have shown that, as we 
see from Figs. 6 and 7, there are such LEC's for which 
the effect of the short-range interaction on the half off- 
shell T-matrix (pi|T(_Ep -I- iO)|p) grows rapidly as the 
kinetic energy of outgoing protons decreases in the re- 
gion needed for KVI experiment. Note that just for such 
LEC's the short-range interaction that is described by 
Eq. (|49|) cannot be associated with any potential. Such 
an effect of including the short-range interaction might 
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seem very surprising, but the difference between the pre- 
dictions of the potential models and experiment in this 
kinematic region grows in the same way. This lead us to 
the conclusion that the inability of the potential models 
to incorporate the short-distance effects which manifest 
themselves when we consider the problem from the point 
of view of the more general approach based on the use 
of the GDE, may be the source of the over prediction of 
the potential models in describing the pp^ cross section. 
As we see from Figs. 4,5,6, and 7, these effects can im- 
prove the off-shell behavior of the 27V T-matrix and the 
same time keep the phase shift unchanged. This allow 



one to incorporate the above effects into the machinery 
based on the use of the high precision potentials keeping 
the phase shifts predicted by these potentials unchanged. 
As we noted, in this case the operator Bshiz) should de- 
scribe corrections to these short-range components. One 
may hope that such an improvement of the description 
of the NN interaction will make it possible to remove 
the discrepancy between theory and experiment in de- 
scribing the pp bremsstrahlung. Taking into account the 
above effects of the short-range interaction on the off- 
shell properties of the 2A^ T-matrix may also open new 
possibilities for explaining the Ay puzzle. 
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